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l. Introduction 


1.1. Purpose and Problem 


purpose ol th ! tigation is to determine 
me oand eo pai cles lound tina variety 
rere al nha experi ental elastomet! latices 
to relate the findin with the ise of such 
or the mip nation of leather as a substitute 
ition by mea of polymer solutions |1 

neans OF WW SILU POLS ‘ ation |2 
: comparison of tl behavior of different 
nth mMmprei on of leather and for a com 
) thre ! pat ( (| stributior _ I Wiis 
hat a nel eproducible method common 
f them would be nee irv in order to discount 
and artifacts inherent in anv one method 
0 t wo i be rece arv to evaluate the 
‘ nd lin ) he method selected 
! ior deter rit ! Lisa pe SIZt nha ire 
{ rit ol ! vinthe ( clastomet! 
estu I \ mae by he methods 
0 diffra ) ittering iriace ad 
' ractionatior nd microscopy Because 
cessitv for data on the shape and mode of 
on, as well a the particle-size distributions 
pure synthetic elastomer s ISpensions micro 
was chosen Inasmuch as the size of material 
he colloidal range the eleetron microscope 
eded to provide the direct means for observa 
ivallable methods of specimen preparation 


on the accuracy ol electron 
diff 
tv of such methods lies in the nature of the material 
ng examined With 
threefold | To 
as individual particles 


pose severe limitations 


roscopical measurements An inherent 
latices, 
maintain the 
material 2 To 

the nonrigid particles in their original, presumably 
pherical, shape without changing their size 6 To 
parat« the 


synthetic elastome 
ihe problem Is 


spersed 


polymer particles from other material 
the suspension 
If one applies the usual method of specimen 
eparation involving direct dilution with water and 
ryving onto a_ collodion the particles 
latten and merge into a hazy layer in which one 
article cannot be distinguished from another. With 


substrate, 


P ; ao 6 et} loft ner 
woke ‘ stu fere ‘ pape 
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i mat i i 
" ul i 600 SOO a 7) 
further dilution, lavers of elastomer of smallet urea 
are formed but still no parti ate structure is 
exhib ted 
1.2. Evaluation of Existing Methods 
Two approaches have been used to overcome ome 


ol these lithe iltie . In one thre latex Was lisp rsed 
1h oath ACGUeOUs solution of polyvinyl alcoho! and 
pon drying, the particle became embedded in the 
resulting plastr film ) lt was found that the 
tendeneyv of the particles to flatten und coalesce was 
reduced lLlowever. this technique Lenppaa image 
sharpne ss and contrast because of additional seatter 
Ing and reduced transmission of the electron beam 
Shadoweasting such a specimen is not effeetive and 
therefore observations of the third dimension and 
the extent of flattening of the particles are precluded 
Also, anv extraneous matter such oa emulsifving 
agent, catalyst, inhibitor, et in the original latex 
ISpension W ll contaminate and partially obscure 


im the 


approach thre 


the elastomer particle : 
In the othe 


brominated |4 


preparation 


lntex Ispension Was 


This served to diminish the amount 


of flattening by hardenin the soft particles \p 
plication of bromine by itself may fix the particles 
us sph res, but severe flocculation prevents accurate 
particl SIZ counts The Ist ol : commercial 


Acrosol OT [6], as a 
in addition to bromination has been 
Although thre 
flocculation, the required encases 
and obscures the individual elastomer particles To 


wetting agent, Daxad—-11 [5] o1 
dispersing avent 
deseribed 


dispe Psilig agent reduces 


excessive amount 
i thorough washing of 
the final specimen preparation is necessary 

Neither the exact washing te hnique nor the effect 
of washing on the particl SIZ 
Regarding the shape of the parti- 
cles assumptions have been made that the partie les 
are either perfect spheres 


remove the « Xcess dispe rsant 


distribution have 


bee li dese ribs d 


1] or ellipsoids [5] resting 
The effect 


has been 


entirely on top of the sagging substrate 

of bromination on thy of the particl 
estimated to be about a 9 percent increase in diam 
eter {4] A 9.6-percent increase 
indirect measurement, has 
No correction 


SIZE 


calculated from an 
reported 7| 

for partial 
This could be 
in the case of particle diam 
eters nearly equal to the thickness of the substrate 


also bee nh 


has been made, however 


embedding in the softened substrate 


a considerable factor 





2. Materials Used 


Th dentifieation ind compo the four 
latice nvestivated are contames I ove tabu 
lat on is. & 10 

3. Specimen Preparation 
3.1. Exploratory Methods 

\ ous methods were u ed ttempts to 
eliminate swelling the pout ‘ ly bromimation and 
using a dispersin agent ‘ ! I their original 
harper ana e and retamin ‘ rete eparations 
between partich \ freeze-di thod was tried 
in which a drop ot a dilut eo pension of 
le eX Was placed upon a oOllodion 0 rere pecimen 
creen, supported by a gla microscope slide The 
if slide was immediately placed on a brass block 
previousl chilled to about min loom ¢ freezing 
the drop. The entire block of frozen specimens was 
placed in the chamber of at RCA LV pe EAL | 
vacuum unit where the ie iblimed during evacu 
ation. leaving thr particles O1 the collodion film 

Another method was tried with latices A and B, in 
which the pH of the diluent was n itched by addition 
of ammonia to that of the late suspension and the 
diluted suspension was aries onto the collodion 
substrate 

In third approach, diluted suspensions were 


spraved by means of an ordinary nasal atomizer onto 


collodion-covered specimen screens Three or four 
shots’’ of the diluted latex were applied to each 
screen from a distance of about mn 

Each of these methods of specimen preparation 


ound particles with 
with all 
that was equally 


well-dispersed 
but not 
reproduce ible 


frequently rnve 
ol the 
Because one 
rood for all the 
to the problem of 
aside even though 
promise of producing n good specimen lor parti le- 
size studies with the electron microscope 


consistently 
i 


SOT latices 
metho 
latices was desired, these approaches 
preparation 


method tried 


were set 


sper men 
showed 


3.2. Method Used 


, 
i 


In order to arrive at a single, reproducible method 
suitable for a statisti al determination of the parti le 
size distributions the four synthetic elastomer 


latices it Was tl erelore, de ided to apply PT | combina- 





tion of various methods already reported Beer 
such a large quantity of dispersing agent is requ 
to prevent flocculation of the particles and 

a thorough washing of the dried specimer 
necessary in order to expose the particles as 
vidual units, concern was felt about the prefere! 
removal of certaimm size ranges ( onsequently 
necessary to devise st veral Speci washing ‘ 
niques and Lo evaluate hot only thr comparat 
effects ol thre washing technique son thr particle 


byt nherent limitations 
for 


part 


| 
iii 


{ 


thr 
Tike 
lntex 


thr 


distribution also 


electron ) i aunhs determini 


mic . 
and sha 
\ sample 


oOscopyV a 


| thre —f 
each of 


pe 
ol 
nqueot 


Daxad 


approx 


{ ; 


i 


with ! 


agent 
ol 
drops ol 


mint I 


this 
watch vlass 
} 


bottl 


percen 


ited s iIspension 


nha exposed to brom 


Irom a of hlauid brom 


i iIspension Was ed 
r 
7 


brominates 
distilled 
about () cn) 


to i rubbe concentrati 


W ‘ r 


I 
percent 


1 
| 


FT| 
Dro 


were plac dl upon ¢ ollodion-covered 


ps ol the final suspen 


= | Hnen seres 


and ai 
time required for the « 
| 


j | 
pita ha 


class 
than the 
In 


support a by 


| 


i 


micros Ope 


n hour onyvel i! 


! 
hum a 


following 


itioned kept 
relati 


he 
cl 


an air-con room 


50 pe reent 


to evaporate 


”) . 


vi 


about (’ and 
Immediately thereafter 
ing techniques was applied 
mens of each sampl of latex 

performed either with distilled wa 


each of wa 


| ] il 


Phe 


I 


oO ¢ separate 


1 
Wasillh 


t 


nt about Si}y 


or with steam at atmospheric pressure 


Jet Washing 


a 


stream of hot water from a wash bottk 
the prepared specimen secre 
The slides were held, 

of about 45° from the vertical, and 


was applied at a near-grazing ang 


\ fine 


laved ove! on 
| 


glass slides specimens 
at an angle 
jet of watel 


This washing was repeated 3 or 4 times for perio 


i} 
) I 


15-see duration 


of about 


b. 


Immersion Washing 


The prepared specimen screens supported by 
glass microscope slide were immersed edgewise 
15 i dish of hot water 
Successive immersions Ww 


{ 
i 


an angle of about into 
about 15 
drawals were repeated 3 or 4 times in the same « 


of water, | 


sec ana 
i 
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c. Pipette Washing 


‘ placed upon the prepared 
alter contact 
This process Was 


Jrops ol hot water wel 
men scecrecns and re moved Ld Sere 
means of a medicine droppet 


ediately repeated > or 4 times with fresh drops 

ter 
i. Steam Washing 

The prepared sereens on their supporting glass 

were placed at an angle of about 15° from the 

eal in a wire basket in the upper part of an 

ed glass column about 40 em long Steam 

passed continuously through the column and 

| the specimens for a period of about 15 min 


e. No Washing 


\ was 


from latex sample 


()) roup ol specimens 
ed with Daxad nna fixed with bromine in 
ner deseribed abo except that a slightly 
concentration of bber was used Even 
I » Washing technique was applied un ae 
preparation vas obtained hie | \ 
h this unwashed . paration Was not repro 
the micro raph can ty used lor comparison 
he other washing procedures used in the speci 
oO fo iter nple A 
3.3. Shadowcasting 
ecImM we | and then dried over 
hose pecimet elected for shadoweastin 
( ited to a p of about 10 mm ol 
y or less Pure ¢ olytic chromium metal 
ere lowly evapo ited from a conical tungsten 


i t placed wu cl manner as to creat 
! shadows vho I ths were } nh ’ 
hye titude of the particle causing them 


4. Interpretation of Electron Micrographs 
uUppearanes and shape 
when the latex \ Is 
v diluted with water and allowed to dry out on 
This particular field offers 


llustrates thre 


Figure 1, B 


particles sample 


he latex 


ollodion substrate 


nformation than the usual micrograph ob 
irom such a pre paration The rubber late X 
enerally coagulated into large flattened masses 
howing no particulate detail at all Figure 1, C 


f the same latex in whi h the diluted 
treatment by exposure 
Instead ol formless 


example 0 
was given additional 
and allowed to ary 


wipe 


rrommine 


masses, round particles are observed. However, the 
wor dispersion of these discrete particles prevents 
measure of their size and numbet 


i irnte 
| vrure 2 
iple B prepared as described in se 
ire 2, A, the specimen was given a jet washing 


represents electron micrographs of latex 


tion 3.2 Kor 


2, B, immersion; 2, C, pipette; and 2, D, steam 
he presence of resid ial mater al in the background 
f these micrographs indicates that the specimens 
were not overwashed, and that only a minimum 
number of particles was removed Excessive wash 


will remove not only the extraneous material but 


ny 

also rT much vrreatel hn imber ol particles The 
delicacy of the steam washing is shown in figures 
1. D. and 2. D. where the covering of dispersing 


avent has not been completely re moved The CXIst 


ence of such an envelope around the particles is born 


out by their diffuse outline and by the shortened 
lengths of the shadows, which would have been five 
times the diameter for unembedded hard spheres 


The envelop ol dispe rsant prov ides an advantage in 
of the surface of the 
particles from coagu 
addition 
Originally 


that it reduces the free energy 
watel helping to prevent thre 
lating during the drying-out process and, in 
cements them to the 


that 


collodion substrate 


it was beheved steam as a method of washing 
would prove to be the iy 31 


tha dispersant leaving 


in that it would dissolve 


all the 


parti les be hind 


uwiy { 

In practice however, in orde! to accomplish adequate 

washing, it was found necessary to allow the steam 
so that the actual washing Was done iL 


Lo conde lis 
Watel 


Attempt 


specimen 


by hot 


drople ts on the 
additional 


and im part 


part by vapor 


to prevent condensation of 


by 


preparation employ ne insulation 
infrared radiation for 


around thr lass column and 
the time required for effective washing destroved the 
collodion film 

The electron micrographs seen in f ( epresen 
the results obtained when the four different laties 
wer ven thre Sai washing technique that 
pipette Th conspicuou differences mn the number 
and the sizes of particles should not be taken to 
mean that these differences result from preferential 
re moval ol parti les In thre wash hye proces Rather 
thev are most likely due to variations in the original 
concentrations pli surface charges, rate of drying, 
the ratio of the substrate area to the volume of the 
drop ol latex Suspension and othe acto that 
influence the settling out of the particles onto the 
collodion substrate No single selected view 
necessarily tvpical of the over-all observations, and 
a statistical analysis ts required before an conelu 
sions can be mia with respect to partich IZe 
distributions 

Prior to determinations of the partic! s1zZe¢ Lol 
these latices, a consideration of the actual shape hia 
to be made The large sphere seen in figure 4 
presumably a mass of molten chromium that fell 
on the dried collodion substrate and subsequently 

what a rigia sphere lving 


shadowed demonstrates 


on top of the film would look like when shadoweast 


in such a way that the k ngth of the shadow ts four 
times that of the height of the partiel Actual 
measurement gives a value of 9000 A for the trans 
verse diameter of the “sphere 

The transverse diameter ts that distance measured 
across the particle and normal to the direction of 
shadowing This value is about 5 percent less 
than the height determined from the reometryv ol 
the shadoweasting technique This discrepancy can 
can result from sagging of the supporting substrate 
under load and a very slight embedding of the 
particle in the film, as shown schematically in 


figure 5. Furthermore, the inherent approximations 
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Ficure 3. One washing technique applied to four latices. 


A, sample A (butadiene-styrene), pipette washed, shadowed 5:1, B, sample B (butadiene-styrene), pipette washed, shadowec 
C, sample C (butadiene-aerylonitrile), pipette washed. D, sample D (butadiene-styrene-acrylonitrile), pipette washed, shadowed 4 
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“a A STEER eh 





—— el ee Banh tothe meal 


ana made in required lor 


If this 


assumptions thie setup 


shadowceasting also contribute to this error 


kind of precision could nevertheless bye attammed 
on measurements made on the small particles 
determination of the diameters derived from. the 
lengths of thy shadows would = be perfectly valid 
and well within the experimental errors of the 
method Such mMeastirement howevel can bye 
grossly in error because during the removal of the 
matrix ol extraneous material Surrounding the 
partic s the hot water softens the collodion sub 
strate not only causing thr film to sa between the 
wires of the 200-mesh supporting screen but also 
allowing the part cles to. sini nm ana partially 
embed themselves into the elatinou collodion 
membrane In addition, as was seen in the case of 
team washing some of the matrix of extraneous 
material remains, thereby partially obscuring the 
particl \ residual film of dispersing agent and 
other materia n the latex su pensior remains to 
form a laver of varving thicknesses on all the prep 
arations This adsorbed laver of add onal matertal 
can be seen not only on the collodion substrate but 
on the latex particles as well. In figure 4 compare 
the smooth profile of the laree sphere produced bv 


a molten mass of chromium, which never was in 
the latex suspension, with the rougher surfaces of 
the small latex particles to which a laver of extra 


neous material dlispe rsed in the latex Suspension has 
been added during the 

There is flattening of the particles them 
of the latex 
vapors does not compl tely fix the 


drving-out 


process 
also a t 
selves, becauss exposure ly suspension to 
particle 
that 


undistur bed 


bromine 
is assumed 
the particle is a true sphere while in thi 
elastomet Suspension This change in shape results 
from flow of polymer during the drying-out process 
the and 
nonrigid particle 


direction 


in its original spherical form It 


when the forces of surface tension of water 
upon the 
im the 
the vertical and ine reasing the horizontal 
The resultant shape of the particle 
flattened ovoid form indicated in figure 5 
made of the 
the thickness of 
laver of extraneous material 


interfacial tensions act 


These forces all act same reducing 
“diameter 

Is probably the 
No direct 
distortion of the 
the adhering 


amount of sag 


measurements were 
spherical shape, of 
or ot the 
in the collodion substrat« 

All of these 
of the 
calculations to determing 
latex particles based upon 
between length and altitude can 
On the other hand the 
portion of the shadow can be used as a criterion for 
measuring the diameter of these particles only if the 
For the 
that 


a foreshortening 
Thus, 


diameter of these 


factors contribute to 
lengths on the 
the true 


shadow substrate 


reometriu al re lationships 
lead to 
widest 


shadow 


serious errors 


metal radiates in effect from a point source 
the Lungsten basket 
crucible and the chromium 
atoms is relatively large. The “polar cap’’ of metal 
built up during the deposition of the chromium goes 
bevond the equatorial diameter of the particle 
The thickness of the “polar cap mav be reduced 
by the with higher 
number lectron 


distances used conical 


serves as a source for 


use of suitable metals atomic 


produc ing adequate 


one leeeia 
cattering 
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the case of the 
this is 


thinner 
9000-A “sphere 

error The effect of the 
tively small and large particles can be compare 
figure 4 Moreover, comparison of the 
the shadows caused by the small particles se 
figures 2 with that of the 
figure 4 that the triar 
shadows associated with the small particles ar: 


with 
larger 


deposits In 


not as 


SCTLO 
“polar caps on 


shape 


and 3 large partic 


makes it clear almost 


characteristic of contact 
a flat surface 

These cumulative errors preclude the measure! 
determi 
Although 


perfectly ref 


an exposed sphere in 


of shadows as a precise means for 


these 
ten hnique dor s not 


vaarticle diameters of latices 
| 


shadoweasting 


the true shape and size of these elastomer parte 
it does, nevertheless provicd a three-dimensi 
aspect to the electron micrographs and ere 


aid the 
diameters of the 


sufficient observer In measu 


the 


ceontrast to 
transverse (horizontal 

Krom the 
of the partir le upon 


tortion from the spherical shape during the dry 


part 
above considerations re varding SW 
exposure lo bromine ana 
out process in the specimen preparation, a corre 
of the observed diame ters should be applied B 
the uppal 


Contamination-growth of 


of these lactors operate to make 


diameters too larg 
particles by electron-beam 


bombardment was 


to a minimum by the avoidance of long expos 


In the opposit SeTis¢ because of the vacuum 
ditions required for shadoweasting and the a 
examination of the specimen In the electron mi 


scope, the particles may undergo a shrinkage thro 
{ the indetermu 
quantitative 


monomer In view of 
nature of these 


could bye applied to the observed data. 


loss of 


errors no corres 


5. Measurement of Particles 


At least 10 micrographs of representative field 
each of the samples for each washing technique w 
taken at 10,000 diameters. The RCA type EMI 
50 kv.) electron microscope was used for this wo 
Calibration of the instrument was obtained by m« 

Formvar (polyvinyl formal) r 
speculum metal diffraction gra 
rulings per inch. This metho 
reliable within 2. percent ‘T 
original negatives of the electron micrographs vi 
projected on a provide further lu 
enlargement of about 30 times \ 10-cm s¢ ali 
0.5-mm graduations on a glass plate was simulta! 
ously projected with the eleetron micrograph \t 
magnifications used, the smallest scale division re] 
sented 500 A. A pair of dividers was used to comp 
the diameters of the particles with the superimpo 


) 


of a shadoweast 
of a planar 4-in 
having 20,000 


calibration was 


screen to 


scale. Some of the extraneous material appeares 
the final preparation as discrete particles less tl 
100 A in size In order to avoid confusion v 
elastomer particles, an arbitrary lower limit was 
( ‘onsequently , no particles under 100 A were count 
The diameters were measured to the nearest LO 
The values reported here are the apparent diame 
of the particles as observed_in the projected elect 
micrographs 
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6. Statistical Data and Their Physical according to size into six groups. The first group 


Significance consisted of the smallest particles comprising 10 per- 

cent of the total count; the second group of the next 

\ comparison ol nveravt particl diameters larger le percent the third and fourth groups of 
25 percent each; the fifth group of the next larger 


tained trom individual micrographs of unshadowed ~ 
i percent; and the sixth yroup consisted of the 


shadowed preparations showed no. significant 
largest particles making up the final 10 percent of the 


ferences No distinction was made therefore, pa . 
ding this particular feature of the specimen total count Ihe range ol particle sizes covered by 
reparation A statistical analvsis 11) was made each ol these percentil yroups IS plotted versus the 
mi othe Viewpoint ol determining what effects. if corresponding group number The median value 
the washing techniques have on the partick occurs at the boundary of the third and fourth pel 
distribution The results are indicated in centile groups Following the usual procedure, it 
The Irequency distribution of particle was considered advisable to omit from figure 6 the 
smallest and highest peres ritile groups because of the 


correspond! a en as co “ | S : 
p ding “ ach wa hing technique considerable variability that exists at the extremes 


imarized for each sample in the following manner. | fo, samples of moderate size. The distribution cor- 


particles of each size, as observed in the various responding to the unwashed specimen of sample A 
| 


rographs corresponding to a particular sample is plotted in the same manner in order to provide a 
a vashing technique were counted and arranged comparison with thre washed specimens ol sample \ 
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The graphical presentation of thes ita shows no 
consistent svstematic differences amo he four 
washing techniques with regard to particle-size dis 
tributiol \ more detailed statistical analvsis does 
uncover some significant vamratio among the dis 
tributions corresponding to any one sample, but the 
lack of consistency of these effects fron sample to 
sample and the smallness of their va tions renders 
them unimportant from a practical point of view 
For sampl \, the effect of washing can be evaluated 
“ th relerence to the imwashed samp | rure © 
shows that th effect revgardiess oO the washing 
technique emploved, is slight 

lable 1 lists the essential statistical information 
relating to the particle size measurements described 
in this papel The medians and averages are based 
on measurements obtained from all the elastomet 
particles counted i he average values obtamed by 
the various techniques for any given latex do not 
differ from each other by more thar percent for 
sample B, 9 percent for sample C, and 6 percent for 
sample LD In the case of sample A, the discrep 
ancies are somewhat large approximately 20 pet 
cent) because of errors of an accidental nature that 
occurred during the preparation ) ing the steam 


washing of this sample the specimens were 


prema- 


turelvy removed from the steam bath and insuffi ently 
washed, making the measurement of the particles 
difficult because of uncertainties regarding § their 
edges Judged on the basis of their standard 


by the swmbol s-. the 
listed show that anv of the fou washing 


could be 


deviations, denoted averages 
techniques 
used to distinguish latices A and B from 
latex C or D, and latex C from D 
It is to that 
median particle size is smallet the corre- 
This ndicate an 
apparent skewness in the particle-size distribution 
Very probably this effect was brought about by the 
fact that no parti les smaller than 100 A in diameter 


in almost all cases, 
than 


would 


interesting note 
the 


sponding average value 
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were recorded. Therefore, the distributions ai 
reality, truncated at the lower end, rathe 
skewed 

An indication of the reproducibility ol the ex] 
mental procedure may be obtained from the m 


parti | 
For t 
computed These 


the var 


fluctuations of 


tude of the 
culated from individual micrographs 


the values s 


average 


post were 
ations character! abl 
the Irom t 
corresponding to a given latex 
nique Thus, the standard 
the fluctuations of avera 
whereas the standard deviations 


tuations of averages based on all 


standard dey zing 


obtained he singel 


averages 
and washing 
deviat 


res of single 


microg 


( xpre S the 


particles frot 


combined total of a// micrographs Accordu 
every 8s» Should be larger than the correspond 
standard deviation s,; by a factor depending or 
ratio of the number of particles in a single m 
graph to that of all mi rographs combined \ 
tistical examination based on an analysis of 
Variances shows howevel that in each case 


the two standard deviation 
larger than can be accounted for by this ratio 

that variability exists 
graphs than be inferred from 
range of particle sizes within a 
These findings in line with the 


made in section 4 about the conspicuous ditt re 


differs nce be tween 


more 
would 


among 


the 


suggests 
nuvi 
nre observa 
observed among the individual micrographs of 
2 and 3 Consequently, the advantage gained f 
the examination of a greater number of microgri 
is larger than that due to the mere increase in 
number of particles. From the 
parison of washing techniques, it 
the additional variability among 
approximately the same for all washing techniq 


viewpoint of a « 
is significant 
micrograph 


microgra 


single microgra] 





| 
{ 
i] 


0 





t} 


On the basis of the foregoing analysis of the da 


it appears that each of the four washing techni 


is equally suitable and there is little reason 














inv sing as being the most 
fo all the latices studied here These 
) ined particles inging in size from ap 
{ 100 to 2100 A The use of anv of the 
echniques resulted u the same classifica 
he 4 latices into separate groups, with 
. ticle diamete! of about 600. SOO. and 
orresponain ot copolymers of buta 
( Dutadiel ) lonitrile and buta 
ne-aervionitrile espectively Alls the 
hniques vielded results that agreed within 
intel » percent ind had auppronr mately 
orce Ol eproducipiit Consequently 
on ot sil Vasi hu qu might ay 
) { o Inasm ich is the 
a ‘ Ss expe liv easier to cari 
ontrol, tl | mia ntativel 
adas prete I 
of the fact that ibsolute valu s lor the 
S and ftreque istribution of these syn- 
tices are not available for comparison, the 
ecuracy of the sults obtained by electron 
py cannot be certamed These results 
provide data avatlable at 
ad sure ! (| it) Intercomparl 
latices 
t ot con that one can draw trom 
s that althoug! data do not permit 
wt quantitative recommendations 
heles I tial qualitative and 
antitat ation regarding the ef 
mus tect specime prepara 
tl listribut particle sizes for the 
SCODL tion of the fo v1 
~ Witl dge of the pore-siz 
l ( the data shown In 
ruil il I th technologist to 
mil es of a given latex 
’ Lt I pore st ot 
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this opportunity to express 
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Second-Order Transitions of Rubbers at High Pressures 


Charles E. Weir 


l. Introduction 


expertmental pressure 
PVT) data were presented for 


previous 


report 


temperature 


sulfur vuleanizates in the ranges 10 to 2S 
nt of sulfur. 10° to SO° C. and 1.000 to 10.000 
No interpretations of these data were made In 


ort In this study the significance of the 


ith respect to unresolve 


i questions concerning 


ond-order, or gla ransition in rubbers are 


(i here iis # 


mite! 


lass transition (denote detines a 


ature or temperat valin which marked 
ves In properties oce in high polymers, glasses 
classy State 
whether the 


representative of a true 


ther materials that may manifest a 


Col siderable doubt ns to 


henomena observed at 
odynami 


ess in which attamment of equilibrium may 


transition or merely of a very slow 


be considered temporally impractical [2 Resolu 
| f these problen s is of great theoretical im 
particularly with respect to the third law 
f thermodynamics but Is beyond the scope ol 
this report The widespread interest in the glass 
transition comes from its practical importance in 
mercial applications, and there is considerable 
literature on the subject such as the extensive 
r vs by Sover and Spence } Berge! 5 
\lorey [6], Kauzmann [7], and others. 
There is no established tsothermal volume change 


that takes place at ro ana 
though ¢ 


similarly no latent heat 
and presumably C,) change markedly 
lo a first approximation it would appeal 
that 7, should be ind pendent ol Th 
interest in the effect of P on i. is indicated by the 
numerous recent attempts to ae duce this re lationship 


from data given by Scott [8 


the retore 


pressure 


Scott [8] previously reported P V7 data on rubber 

lfur vulcanizates over essentially the same tempera 
ture interval but to a maximum pressure of only S00 
bars These have represented the only 
available data on high polymers with a direct bearing 
on the relationship between P, V, and 7, 


at the time the measurements were made the future 


heretofore 
Howe Vel 


significance of the discontinuities in slopes of isobars 


vas not apparent, and the data are given only by 





means ol graphs ana empirical equations that are 
ult to use With subsequent increas 
ing interest in the nature of the processes r¢ sponsible 
for the changes observed at pf the data of Scott 
attention trom recent 
view of the way the 


somewhat diffi 


have recelved major workers 
It would appear, in 
presented and the fact that ther future significance 
anticipated 
have derived 
unsoundness of such cone lusions Is con 
firmed im instances by the origmal numerical 
experimental data that have been 
through the cooperation of A. H 
previously high 
essentially the 


Scott s 


results are 


that several que stionable 
recently these 


could not bye 


conclusions been from 
data The 
most 
available 
The 
pressures [1] are 
sulfur and 


and to give 


rrnancle 
Scott 4 
reported data at 
designed to 


Covel SiLirie 


temperature ranges studied by 
wn UNeEqulvor al answer to the interde pendence of 7, 
and P matters of interest her 
Rubber-sulfur that 
several specimens were expected to pass through / 
in the interval although 
were made of specimens exhibiting no 7, in this ten 
These were added to round out 
in SO far as possible 
could not bye 


as well as othe 


vuleanizates were selected so 


temperature ised studies 
perature interval 
the pieture 
sulfur 


bye Cnuise 


However very lov 
studied 
action of the confining liquid a ligh 


The nu 


empiri al equations 


eontents conveniently 
of the 
petroleum distillate on 
merical data obtained, and the 
derived to fit the data, have been piven pre viously l 

it must be ¢ mphasized that subseque nt interpreta 
tions of these data can be considered to apply only to 
the pressure interval studied (1,000 to 10,000 atm) 
and to the rubber-sulfur svstem of high 
sulfur content 


such specimens 


vulcanizate 


2. Results and Discussion 
2.1. Isobaric Transition, T, 


calculated from the numerical data are 
shown for six vulcanizates m figures 1 to 6, In which 
specific volumes are plotted as ordinates. Each 
isobar is labelled with the corresponding pressure in 
units of thousands of atmospheres. In these figures 
the isobars are drawn as linear above and below / 

which is taken to correspond to the intersection of 
these lines. In figures 1 and 2 no transition is ol 

served, and they are included to illustrate other mat 

ters to be pointed out later In figures 4 and 6, the 
well defined: in the 


lsobars 


transition 


temperature is not 








former because of scatter ¢« ffpoints and in the latter 
because only point occurs at temperatures 
above 7, In figures 3 and 5, 7, is considered to be 
sufficiently well defined to concluds 
sideration of the experimental errors [1 


a single 


with due con- 
that there 
is no evidence for any measureabk in 7, with 
It is suggested, therefore, that 
isnota 


change 
changing pressure 
for these materials in this pressure range, 7, 


function of pressure to a rathet rood degree of 
approximation 

It will be noted that some discrepancy exists 
between the values of 7, given here and those of 
Scott ‘ from which these data were essentially 
derived These discrepancies are to be attributed 


to uncertainties in deriving the specific volumes from 
the graphical data of Scott |S} and inherent varia- 
ations to be expected in lime s drawn through experi- 
Therefore mav be 
values shown here. 


significal 
of T, 


to recognize that ho 


mental points no ce 
attached to the 
It is of the utmost 


factor appears to act on compression 


actual 
Importance 
to produce any 
measurable change inh / 

Specimens used in these studies will ultimately be 
analyzed and studied dilatometrically atm No 
significant change in any conclusion presented here 


x 


of 
a 


is to be exper ted as a result of anv measurements at 


latm. Slight revisions of specific volumes or values 
of 7, 
due course 

From the fact that the isobars represent a family 
of divergent lines, shown best in figures 1 and 2, it 


possible to deduce qualitatively the effec ts of pres- 


mav be expected and these will be reported in 


Is 


sure on expansivity and of pressure and temperature 
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TEMPERATURE a 

FiGuRE 1 Isobars for rubber-sulfur vulcanizate containing 10 


percent of sulfur 
The number opposite each isobar denotes the pressure in units of thousands of 
atmospheres 


312 


Such effects 


on compressibility interrelat 
the geometry of the family of lines 
the fact that the expansivity is proportional 
slope of the isobars and the isothermal « ompress 
to the separation of the tsobars at a wriven ter 
ature Qualitatively the expansivities 
pressibilities behave in a 
pointed out by Bridgman for solid 

From the empirical equations derived previo 
the compressibility may be written as 


are 


isobar 


ana 


manner similar t 


S and liquids 


1 OV 
-la—bP+cP a’—b’P+e'P 
V,oP | 
a’’ 5’ P c’P 
where the coefficients are essentially positive ly 


empirical equations it is noted that 6’ and ec’ 
change signs rarely; but a’’, 6’’, and e”’ are foun 
reverse in sign somewhat more frequently 
change in compressibility with temperature is th 


1 OV 


b’ P+. 
V, 0Po7 \ f 


‘'P?+2(—a’"’+6”"P 


c’P*t 


with similar considerations regarding the signs of t 


coefficients. At a given pressure it Is apparent fr 
1) that the compressibility increases with 
creasing temperature, since a’ > a” ,b’ > |b” 
c From eq (2) it appears that at a given pres 
the rate of increase of compressibility with temp 
ture is less at higher temperatures that 
T,) than at lower temperatures (that is, below 


eq 


Is 


abo 


j 
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However, for some specimens exhibiting 7, in this 
temperature interval, a’’ is that the 
reverse is true. This fact, together with the errors 
inherent in defining the smaller coefficients, renders 
variation of compressibility with 


negative so 


a decision on the 
temperature somewhat uncertain. From the equa- 
tions it may be stated with confidence that 
any change in compressibility on passing through 7, 
From 
such an 


however 
is not large figures 5 to 6, 1t can be con- 
cluded that compressibility 
would probably decrease with Increasing pressure and 
that such an effect, if it exists, Is low 
pressures Consideration of the experimental data 
furnished by Scott [9] vields the conclusion that 
existence of a change in compressibility on passing 
through 7, is very doubtful, and if present is certainly 
Jecause compressibility is generally 


Increase ih 


large! at 


rather small 
not determined with the precision of expansivity and 
Is a strong function of pressure and temperature, any 
small such a 
subject to some There is 


( hang Ss in 
uncertamty 


conclusions involving 


property are 


little doubt that a change 1k ompressibility of the 
order of 2 or 3 to 1 Is not shown by these data 
Inasmuch as eq (1) and do not readily permit 


an unequivocal general analysis, it is of interest to 


show the compressibility-temperature behavior nu 


merically Such calculations have been made for 
several vulecanizates with results similar to those 
shown for a flexible rubbe percent of S) and a 


rigid rubber (22 percent of S) which shows the marked 


transition of figure 5 The results of these 
calculations are shown graphically in figures 7 and 8 


isobari 














in which compressibility is plotted against tem- 
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perature, the pressure being denoted on each 
in units of thousands of atmospheres. It w 
noted that a discontinuity in compressibility at 
absent at all figure 8, 
bility isobars being essentially as smooth as 

for the 13-percent vulcanizate shown in figur: 

Differences between the shapes ol 


pressures inh the comp! 


comparison 


curves in these figures are not of interest her 
principal point of interest being the essential smo 
ness of the curves It is to be noted that 


assumptions involved in calculations of the o1 
data 1] have negligible effect on 
which are defined essentially by the shapes of 
isotherms and not by their 


Discontinuities in compressibility 1 


these cu 





compression abso 
magnitudes 
show up as irregularities Ink the original experime! 
data irrespective of any extensive properties of 
material No have 
within errors ol 
measurements 


irregularities ever | 


experimental 


such 


observed the 


It is of interest to calculate the change of /- 
increasing pressure on the assumption that the 
nomenon is representative of a true the rmodyna 
The equations def 


transition of second order 


this variation are 





ip? / ( ('p x y 
and 
i}? if a ¥ J ) 
where the terms have the usual significance 
Equation (3) has been evaluated by Gee [11] ‘ 


Sufficient data 
From eq 


found d7/dP deg/ bat 
available here to evaluate eq } 


re] B. AB 1s found to be 


0.016 
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( and ¢, refer to temperatures slightly below 
slightly above ? Krom the ¢ mpirical equation 


, Ae 1s found to be 


Aa=2(a a’ P—6"’ P?+e’’ P*)(t,-t 6} 
bv substituting q (4) mav be shown to bs 
P—(q’—b’P a’ +b’ P p 
f Ia h’’pP { j? 7 
hn reduces at low Dt! res to 
ridP ! 2a X 
( s1dering the 22 percent rubber sulf ll vulean 
n parti ilar although all other rubbers gre 
comparable f fo TAY bye taken to be appronri 
\ 10°. whilea 3 6 oO a’ iss1.3 10 
s 1.4 LO Che value of d7/dP is then 
! order of 7 () alt atm This corresponds 


change of 7, of the order of 7 dee C/1.000 atm 


that the 
denominator of eq 


yw pressures It orders 


Is not apparant 


tude of numerator and 


o not change with increasing pressure How 
simpli numerical computations using the 
ents of the empirical equations show that 

his the case On the basis of these data, a chang 


. of the order of 70 dee ( is predicted from eq } 


om the value given by Gee 11] for eq } the 
esponding change would be approximately 160 
( 10,000 atm It is obvious that no such change 
\ s in figures 3 to 6; and indeed, any change in 
th pressure il presel t al all would appeal to be 
than | deg © 10,000 atm The discrepancy 
predicted and observed change in 7, with 
¢ /? would appear to lend support to the sug 
ons that the ph nomena occurring at Y are not 
sentative of a true thermodynamic transition 
ond ordet 2 } 
] 2.2. Isothermal Compression 
lypical isotherms are shown in figure 9, and 
ise all other isotherms art essentially similar 
here seems no reason for their reproduction here 
All isotherms are inherently smooth, and it must bs 
oncluded that a second-order, or glass, transition 
ssure Fo does not exist in these data The 
stence ol ie would appear to be of Importance 
om both theoretical and practical viewpoints and 
possible reasons for failure to observe such a phe 
menon must be considered 
First, occurrence of such a transition at low pres 
es, that is, below 1,000 atm may be essentially 
led out, since Scott [8S] observed no such transition 
- mn to S00 bars 
90 


Second, 


the effect of the confining liquid on the 


unsition may be discarded as measurements of the 


mount of such liquid absorbed vielded results ot 
ich small magnitude as to invalidate this possibility 
Third, the 
ompressibility 


change of 
might 


small magnitude of any 


occurring on traversing P, 
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render observation of such a phenomenon most diffi- 
cult, particularly since the compressibility is strongly 
pressure dependent The difficulties from 
this possibility cannot be minimized and are very 
real However, it would appear from the verv real 
similarity of the compression curves of all polymers 
without regard to ea or the che mil al 
structure of the pelymer that it is very improbable 
that any such 7, exists 


arising 


studied to date 


(ll previous studies ol this nature have re ported 
smooth /7| Adams and Gibson [12] 
ised essentially the experimental apparatus 
applied here sridgman [13] | 
different method requiring no confining liquid and 
has reported some irregularities. Tammann and 
Jellinghaus [14] studied glasses at various tempera 
through the 
? OOO kg em 
no isothermal irregularities 


relationships 
Sarre 


iowevel has evolved a 


range al 
re ported 


tures transition pressures as 


high as and likewise data 


showing 

The report of Bridgman deserves special considera 
The data from which this report 
were made available through the cooperation of P. W 
Bridgman [15 Krom these data it would 
that the isotherms themselves are also smooth, but 
irregularities appear mM the differences of the exper 


mental data Jecause all rubbers studied which had 


tion originated 


appeal 
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similar rubber content did not exhibit similar behav- 
irregularities 
lhis is not 
considered 


ior, it appears most probable that the 


did not arise from the polymer alone 
typical of the 


believed to be transitions 


here 

It Is concluded therefore that there is no rood 
evidence for a /’, comparable to 7, No absolute 
requirement for such a 7’, exists because /’ and ] 
are independent and observed behavior in the \-7 


plane does not imply a definite behavior in the 7’-| 


plane Howe ver, sine | Is col side read to bye of vreat 


importance in determining the behavu 


phases and values of \° smaller than those 


rofl condense d 


‘ x per ted 


at 7 () I by thermal contraction alone can be 
produced by these pressures it is of interest to eval 
iate internal energy changes duc to isothermal com 
pression over this pressure range and by isobari 
expansion over a temperature interval that include 


Y,. These 


- energy changes may be evaluated by the 


easily derivable re lationships 


Ak 


Ak 


where the quantities have the usual thermodynamk 


significances These equations wer evaluated for 
natural rubber in which 7, ~200° K he tempera 
ture limits were 300° and 175° K, which include 7, 
and the corresponding pressure limits were 1 and 
10,000 atm, respectively The following assump- 
tions were required in evaluating eq (9 and (10 

1. 8 for natural rubber was assumed to have the 


Same pressure depends nee found for the 
sulfur vulcanizate studied 


1O-per ent 


2. a was assumed to be a linear function of pres- 
sure, with a pressure dependence similar to that 
found for liquids by Bridgman [16] 

3. Quasi-equilibrium was assumed at 7, so that 


10) could be applied 
out by experimental data [17] 
to be 


eq This &aSSUMpPLLON Is borne 


1. a was assumed continuous in a mathe- 
matical sense at 7, 

Using Bekkedahl’s data for specific heats [18] and 
expansivities [19], together with the above assump- 
tions, the following values are obtained 


AF), 103 j/g 


(AE), 200 j/g 
These values are of the same order of magnitude for 
natural rubber where 7, is far below room tempera- 
ture If harder rubbers are considered, the lower 
limit of temperature required on the integrals of eq 
10) to encompass 7, becomes higher, and for very 
hard rubbers the value of (A/)p reverses in sign as 
T, becomes greater than room temperature. The 
value of (AZ), however, will remain relatively un- 
changed. It is clear therefore that for the rubbers 
studied here that |(A/)-7| is certainly not necessarily 
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Ak); Therefore, a quantitative 


smaller than 
ference ip internal energy changes cannot be uss 
account for the absence of ?, on isothermal com 
It must be emphasized that the source of 
respective internal-energy changes has not be« n 


sion 


sidered. On isothermal compression, it is most } 
able that the decrease in internal energy arises | 
cipally from change in potential energy [16] 
in isobaric contraction 
uted principally to kinetic sources. It appears 
sonable to suppose that the transition is not ca 
by potential effects, particularly 
theories of the causes of the transition in poly: 
+, 17] emphasize the retarded viscous-responss 
cept The latter of these, deduces 
viscous-response concept from an hypothesis ol 
ical free volume The free-volume concept we 
appear to require some modification becauss 
volumes attained in the present studies, altho 
considerably less than those expected at O° K 
thermal contraction alone, produce no tsother 
transition and appear to have no effect on the isoba 
transition. In the light of the present evidencs 
extension of such concepts appears justified 
Important quantitative conclusions differing fr 


, whe 
the de rease 18 to be att 


bec nuse rec 


howe ver, 


those of this report have been deduced from the data 


obtained in previous studies at elevated pressur 
Although the measurements were made by Scott 
and Tamman and Jellinghaus [14], pertinent con 
sions of interest here have been reached by la 
workers from a study of the reported data ‘ 
conclusions are frequently erroneously attributed 
the original authors, and, in 
to be highly questionable 


most instances, 
It is necessary, therefo 
to examine those previous data critically and 
evaluate conclusions drawn therefrom tl 
differ from those of this report 

The frequently quoted results of scott [8] W 


such 


studied rubbers to a maximum pressure of 800 bars 
are ¢ ited chiefly with respect to the behavior ol the 
under 
‘'s a glass transition in the temperature interva 


5-percent rubber-sulfur vulcanizate that 
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rted The original experimental data on this nterpretation, show no change in T, the same pres 


] canizate were made available by A. H. Scott [9], sure interval 
| are reproduced in table | From these data the Numerous qualitative and quantitative deductions 
pressibilities AV/VoAP were calculated by have also been drawn from the data of Tammann and 
differences These compressibilities are given Jellinghaus [14], who investigated glasses at 100 
table 2, and a plot of compressibility versus tem- | kg/cm?’ pressure intervals to a maximum pressure of 
ature is shown in figure 10 for interpolated pres- | 2,000 kg/cm Isothermal measurements were made 
es of 100 and 700 bars No appreciable discon- at approximately 10-deg-C intervals cove ring a 
uty in compressibility is manifest at either pres range that meluded 7,. These data have been 
re, and similar results will be obtained at other | widely quoted as showing a pronounced change in 
ressures as is evident by an inspection of the data | compressibility on passing through 7,. The authors 
f table 2. estimated theu experimental error through a meas 
The generally accepted conclusion that Scott’s | urement on water that was compared with the early 
data show an increase in 7, with increasing pressure 1893) data of Amagat [20 By such a comparison 
! s also subject to question The isobars given by they report an average error of 3>10~* ml/g and a 
. Scott [8 are insufficiently detailed to warrant definite maximum error of 1.1 10 ml/g in specific rolume 
onclusions. Jt must be emphasized that if the widely Cursory examination of the specific-volume data 
toled 165 deg C rise in T, at 800 bars for the 19 5-per- does lead to indications of a marked change in 
' cent vulcanizate is literally accepted, one must be | compressibility at 7, as widely accepted. However, 


repared to wnore completely similar data jor the study of the isothermal first differences of the experi- 


ys 


-percent vulcanizate, which by simular rigorous mental data reported discloses a considerable scatter 
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values 


Al 


The first differences to the 


eorre spona 


experimentally measured, and it is these values 
which determine compressibility The estimated 
average error of 3>¢10~* ml/g in specific volume now 
amounts to a large proportion otf the AV value 


These considerations apply to the isothermal data 


and if it is now required to estimate the temperature 
coefficient of compressibility, it Is necessary to com 
pare the Al values for a given pressure at the 
arious temperatures The change in compress! 


bilitv at 7, is precisely this temperature coefficient 
of compressibility that ts clearly subject to consider 
ably more uncertamty tha tl compressibility 
itself. To illustrate the problem, the first differ 
ences of the reported data on salici 3/14], which 
represent Al values for each 100 ke cn mre plotted 
wainst pressure for four temperatures in figure 11 
The vertical dimension of each poimt is taken as 
; 10) mig im accordance with tl thor’s esti 
mated average error | pec ) ind ms ind 
cated by the barred line In a statistical sense this 
estimate for error un the dill ces probably con 
siderably too low, but all perti t is1ons are 
obvious, even with the underes ( ) 

Thre erticau separation of hy peo i viven 
pressure will be a measure of the efl ) mperature 
on compressibilit It s to Ie inticipated = tha 
compressibility. will increas vith wu ising tem 
perature However nh onl S 0 | 20) pressures 
shown is the order of the pou order ol 
increasing temperature The o ises contain 
at least one nstance of an mnversion micating at 
least one unstance of greatet ompressibilit at a 
lower temperature In any event, inspection of the 
amount and degree of overlapping ar nversion in 
figure 11 indicates the questionable nature of any 
conclusions concerning compressib s related to 
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temperature It mav be noted that instance 
inversion of compressibility with temperature 
occur in the data shown in figure Ss The amoul 
such inversions are of a much smaller ordet 
most of those indicated in figure 11 and are w 
the experimental errors of the present measuren 
A tabulation of compression values (A\V/\ 
three intervals (1 to 1.000: 1.000 to 2.000: an 
2 000 ke em is given lm table > Aas caleulated 
the reported specific volumes of the three 
; we yeas w40e ~ 
. 
4 
i] ” LJ : . ‘ x 
mm” { 4/ { { 
PAL ws w4ue 1 








14 These data show that any conclusions 
ding the change in compressibility at 7, must 
ialihed at least as to the pressure range and in 
of the results of figure 11 should be drawn with 
eme caution 

in be concluded solely from a detailed examina- 
of the data of these investigations that the inter- 
tions, which have been subsequently advanced 


mnection with questions concerning the effect 


ressure on the class Lransition, are not inmequly 


established by thes uta It is therefore 
ested that as a result of this and the previous 
s that there is no valid evidence for a discon- 
ol compressibility ut Ra nor for any appre 
effect of pressure J. in the rubber-sulfur 


m of high polymers 
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Thermodynamics of the Rubber-Sulfur System 
at High Pressures 
Charles E. Weir 
‘ j , ’ r 
T ( inn 0.000 atr . a IX , nal . 
i ( i ‘ i | 
i it . NMi-aty i " ~ 
" 2s l " i 
l. Introduction with respect to pressure are believed to be known 
to a Ssatistactory degree of accuras 2 ill tempera 
oO previous rept 3 4 it have been presented tures It is clear, however, that the derivatives o 
/’ measurements o1 ibber-sulfur vuleanizates volume with respect to emperature are spe ficall 
inges 10 to 28 percent of sulfur, 10° to 80° C known only at one pressure, atmospheric pressure 
© 10.000 atm the empirical tsobart At all elevated pressures eXpansivils must b 
il and P\V7 equations that represented thy deduced from compression measurements and are 
and some u DI ons of the data with therefore subject to all the errors included im th 
» the tsobarn lass transition temperature COTMpression measurements Phe eXpansivitlies are 
cL othe sothermal la transition pressure consequently subject to some uncertaint id an 
2 The preset report describes additional | Calculations dependent on the expat mus 
of the exp: , ita [1] in ealeulation | not be considered as numerically much more signifi 
vf mic qual es of nterest These cal than for order of magn ale 
ms are patternes fter those earned out It was hoped that the en pirical PY‘ quation 
} man im mar dies on liquids \ resume derived [1] would suffice for the calculations deseribes 
» obtained | ! ven by Bridgman | here and thereby reduce the calculations to ordinat 
1) 1 are pres ‘ for the following mathematical manipula ons This expectation wa 
sas functions of pres re and tempera rr not realized and ill data to be presented vere 
ngves f o above for thre byher obtained b iphical method The fact tha thie 
: le ates v4 compression: 2 heat of emp rical equatior Tria not be ised in such caleu 
0 chan n mternal energy on con lations does no eflect doubt on the validity of thy 
mad 4, differer between specific heats quations because the vere derive to reproduce 
ress ! ant volume 0 the volume at iio pressure nna tel 
aut : There is no a prio renusonil tha predicts 
| T | ' 
2. General Considerations oe — of the Feet _ 5 eh ne ilso reproduce 
. Kron hie ab lated OLUTNE char ‘ I VOLUN 
wa onan ons are perunent WO the | a¢ 21° ¢ en earlic ind the expansivitic 
: dies ¥ \ measurements carried ae! ed from. the data ot Scott } th pr ic 
were re ead sothermalh ana olume was calculated at eve pressure am temper 
measurements O elative changes of iture measured Krom the tabular data sobal 
various p The vo e was vere drawn to a larg scale at OOO atm 
it one temperat e (2 ( n il volumes ypical sobars of this nature ha aT el 
ther temperat ng derived from. the earlier [2 lhe isobars may be considered as linea 
ed volume at 2? ( ne ndependent expan within the experimental error of the measurement 
ita This proces has been describes n and trom the straight lines drawn through th 
Theoret call ch measurements coi experiment il pomts the sper fie volun Vil ru 
ee all ‘ 5 ae ; lated at 10-ce nte ils 
‘ iermaily a i ous emperatures con Calculations re , here fo; thre< eubhes 
th a measurement ol expansitity at a single containing 10, 18, and 28 percent of sulfur Phi 
ire sufficient to define all states of interes sulfur contents indicated correspond to the amounts 
In practice howe the quantities required added in. the milling operation | Vulear atiolr 
culation are the vatives of the volums was effected at 140° C for 17 li The f of thes 
spect to temperature or pressurt By virtue compounds s soft and the last on s hard at al 
method of measurement, derivatives of volume temperatures studied; the  18-pereent-sulfu if 
canizate exhibits a glass transition at approx mately 
I iS ( For the 10- and 28 percent-s fur vulean 
zates, therefore the linear tsobars are vell define 
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because all experimental points may be consi: 


the isobars For the 18 percent-s 


two intersecting lines were drawn 


in drawing 
vulcanizate 
slopes of these lines are not as well defined be« 
fewer experimental points determine each line 

glass transition temperature, 7’, is taken to | 
temperature at which intersect 

tabular specific volumes in this ma 
are given in tables 1 All subseq 
from the 


these lines 
obtained 
inclusive 


derived directly 


to 3. 
calculated values a 
of these tables. 


re 


The expansivities indicated in table 3 are obs 
to vary with pressure in accordance with the fin 
of Bridgman In tables 1 and 2 


3] on liquids z 
probable 


that the 
SOTHeW hat 


very dependence of 
Pansivits This 
results directly from the experimental data, wl 


taken literally indicate a negat 


pre ssur 


Too lara cone} 


Is 
‘ eXpansivity in 
instances at high pressures at the lower temperat 
This negative expansivity may be real (see Brids 
3}); but since the experimental errors are such ; 
doubt as to its validity, the linear isobars 
a positive expans 


cnuse 
drawn so as to always show 
the value of expansivity thus obtain 

the influence of the points contrib 
the negative CXpansivity In this connection 


Hlowevet 
lowered by 
it} 
be pointed out that the empirical equations pr 
will also \ ield neg 


ously derived 1} gative expansi\ 
in certain ranges of pressure and temperature be« 
they fit the experimental data quite well. This f 
represents the principal reason why equa 
for these 

Typical specific volume-temperature and sp 
volume-pressure relationships are in fig 


l and 2 


were not considered suited calculation 


shown 


Compressibilities to be used were also cal 
directly from the data of tables 1 to 
first. differences The values so obtained appea 
accurate and as two graphical met! 


that were attempted 


3 by meat 


consistent 


All calculations were performed in experin 
em*-atm The final results were then 
using the 


units ol 


conversion relatior 


verted into joules 


1}— 9.869 cm’-atm 
| A} } ‘ f 
f 
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3. Work of Compression 


ork ot compression denned as 
w=] P| ae 
y 

ited DV graphica ntegration of the pres 
me data ot ib to by mechanical 
The values obtained are given in tables 4 
These data rep I ork done on the speci 
the external compressive forces The work 
ved to be only slightly affected by tempera 
strongly press dependent. Graphs of 
x inction of pressure show the approximat: 


res noted by Bridgman [3 


press 


appear [rom these data however that too 
rvature still exists to permit successful uss 
ir relationship for <trapolation The 
of compression decreases with increasing sulfur 


aires thy 
sulfur 


nt as expected Such decrease is dus 


the higher 
ted that mathemati 


rrations of the appropriate partial derivatives 





lower compressib litv of 


ent rubbers It is to | 


empirical equations {1 vield results in good 
ment with the data ol tables {to 6 that were 
ned graphically Thus the empirical equations 


both the 


resent 


compre s 


considered to rep 
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sibility and the temperature dependence of the com- 
pressibil tv fairly closly It will also be noted im table 
5 that no discontinuity in the work of compression 
occurs at ie This result follows directly from the 
previous report {2} m which no discontinuity in 


ow 


was observed a 


4. Heat of Compression 


; f | 


compression detined a 


r| (=~) IP 


ilated | i graph 


O/7 


Wis cal 


ol / O\ mechanical methods 


by 








Phes 
mip} 
il ho 
ibles 
lr pai 
XY paris 
for the 
pr I 





o obtamed are given in tabl 


mens when 


pre | Kpansivities 

wh to a high accurae and the data ol 
to & must not bys nterpretes too literally 
la the apparent excess crease Ol 
vith nereasi ) =s 3 esponsiblk 

vy values of heat shown u il S at low 
Ih tibial il ! is | i also ob 
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n table 8 at temperatures abovs I, s real i < ( 
ollows necessarily from the increased expan 

bove 7’, 
yical data for wor ind heat of compression 


Internal Energy Change on Compression 


nternal ener: chan on compression 
by subtractin corresponding values ol 
compression from work of compression, con 
both values iis nti nsically POSILLVe The 
obtamed are show?! ! tables lO to 12 An\ 
incies noted mre ts obtained from sub 
data given here n tables 4 to 6 and 7 to 9 
data ot tables Oto 1l2 al ses from the fact performed on compre oO} At som , ' Dre 
ll calculations were performed initially in aure. howevs this proce mu eV e because the 
ntal units of I atn before converting forcea must ultimate become repulsive There 
ts reported her Considering first. the Coe be to he expected that a oe Guan wont 
of tables () gre 5 hich may be considered t] = ; , , . rr r) nal ne chanws 
entially numerica correct for reasons dis oad] tt ; ; snd become positiv: This pres 
eariiel t ws noted that on tsothermal com ie all veras wnnot be reliab estimated on the 
! e internal I | aecreas ina ha nn , o be noted. howeve that 
ated pressure starts to mecrease with in the p sure of reve 1] on the basis of these data 
pressure lr} behavior has been notes ncreas th, . . mp " ‘ i puired 
briden n who | off red the : p nation 1) iby Thiel priaagia Or wane iso) nerves 
er ordinary conditions the molecular forces vith inereasine hardness of rubbet The data of 
cdomnacely 1 o that on compression table | ure more nterestin The marked change 
berated mitia 0 rcompensates thi work mn mternal ene change on pa Nn through the 
herm ns J » be notes Th change 
while ! ! ’ rhe I ol Nn ‘ onabl 
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the lower temperatures is highly suspect as men 
tioned previously as the expansivities in this region 


of temperature appeared to possess anomalously high 


There Is no doubt hnowevel 


pressure cle pe ndene 
nternal 


that at temperatures below 7, th energ\ 
change reverses sign at a lower pressure than at 
temperatures above 7,. This result also must 
arise by virtue of the change in (OV/07)>p at 7, 


Typical data for internal energy change are shown 


6, and 7 


in figures 5 


6. Change in Heat Capacity 


The change in heat capacity with pressure m; 
calculated the relationship OC'p Ol 
l(orV 07 Because it was pre viously ass 


Irom 


that (O1°707 was constant within experim 
erro! such calculations have no basis in 

thes experiments It may be concluded 
assurance that any change in C'p with pressu 


small 
ifference between heat capacities at cor 


quits 
The 


pressure and constant volume mav be evaluate¢ 


the relationship ¢ ( P(oVio7 | 
which involves only first derivatives Full da 
this difference are not given here Partial 
for two rubbers of interest are given in tabl 
and 14 in terms of specific heats. In table 
typical behavior ts observed with the differ 
decreasing rapidly with increasing pressur 
approaching a neghgible value at high pres 
Isobaricaliv littl temperature dependence Is 
which follows from the fact that the expans 


2s of temperature 


Was take! 


] mdependent 
table 14, similar 


above Ol 


behavior is noted 
- but on passing through fa isobarically a mn 
change in ¢ (', is noted. This result is a di 
consequence of the experimental finding that 
appreciable change in) compressibility occur 
# that is analogous to the change in expansl 
The change in the difference between specific | 


act 


heat 





7 
| 


{ 





Oh passing through 7, amounts to approximat 
0.5 }g It is probably fortuitous that this 
agrees closely with the change in specifi 
constant pressure [ol natural rubber observes 


Sekkedahl [5] on passing through t. 
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is assumed that ¢ . 
reasing 


as shown 1n th 


hanged 
may bye 


not mate rially 
further 


following 


pressure conclusions 
treatment 
in the difference of 


ises from a change in ¢ 


} bye shown that the change 
heats probably il 


(’, and that C’) is, therefore, very probably 


ed on passing through ?¢ To show this 
vative of the difference of specil heats with 
to temperature Is take! and is found to be 


a)! 


OV 
7)./(ar) 


| | 
(spor) (SP), | 


the tern m the 
ed that in a sma temperatt 


has bye en shown 


temperature coefficient of compressibility 
compared to the compressibility, so that the 
m is also neg! rie It is known, however 
temperature coefficient of expansivity 1s 
this temperature interval so that the ratio of 
o1/o7 of the order of 4 or } in this 
Pherefore any chang ol f ( ol 
through 7 arises predominantivy from a 
h eXpansivity Because the basic distin 
a ( nad ¢ olves this expansivity 
s that ¢ probably unaffected by 
through ( lol hi normal small 
, to tn expect ‘ ‘ ‘ Irom temperal ire 
Lion Altho li ! Sugvested constancy 
passing through must have significance 
is the processes sponsible for the changes 


this behavior 


ho possible explanations o an 
posed at this time 
the data presente s possible ocak inate 
s of hie ! lat onship 


)/ 


cy = (35) =(35) +(S5) (Sr) 


/ s the internal energy All terms on the 
le of this equatiolr ire obtainable from the 
data However, attempts to verify the 

ted constancy ot ¢ nh passing through I. 
ilculate Cy for any specimen wer nsue 


reasons for this are to be 


The 


ted to the 


failure 


predominant effect of the expan- 


Vv integral in determining the values of AZ, and 
expansivity itself in determining the values of 
7 As pointed out previously, the expan 

ty s not believed to be sufficiently well deter 
d to justify reliance on the results of such 
ations 








7. Pressure Coefficient 


cal 


The pressure coefficient OP /o7 y may be 


culated from the ratio of expansivity to compressi 
bility There seems little point in reproducing thes 
values which are generally of the order of 10 to 20 
atm per deg © It is to be noted, however, that by 


the 


Sufl 


change 


of 
coefhicient is <iscontim 


Virtue Ink CXPANsivIty pressure 
ut / 


Lous wient data 
| 


are presented here to permit the calculation of these 


values if they are of interest 
8. Discussion 

The data presented lol thre typical rubbers ol 
high sulfur content are believed to be representative 
of the P\V7 data obtained her Calculations of a 
similar nature have not been carried out on all 
samples tested due to the laborious nature of thre 
calculations and thre fact that the othe samples 


vielded quilt simular experimental results The 
calculations deseribed above have been carried 
through with one other rubber containing 22 percent 
of sulfiu th results similar u all important respects 
to those obtained on the 18-percent-sulfur vuleani 
Zale deseribed here Specit qualifications con 
eerning the numerical accuracy of the data have 
been emphas Aue but significant behavior of thy 
calculated quantities has been shown to be a direct 


hoOoms W he rt 


ilts 


consequence mental observat 


Is Involved 


obtamed Bridgman for solids and Hatids 
While the general trends of most quantities reported 
are believed to be of interest would appeal threat 
the most significant result lies in the suggestion that 
( very probably chon hot change discontinuously 
il ; Inasmuch as vnificant differenes may 
eXIst il lass-tral Liol cmipera ound hy 
Scott } unal an thee present data t . planned lo 
measure the Dp I Olu! temp i ‘ lation 
ships of the sp nens n these studi ! 
dilatomet« i itl ind modifications of the 
present high p oo pment at el ated pre 
sures The values so obtained will be ised to recal 
culate portions of the data presented here Any 
significant differences will be reported at al re 
ciate 
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l. Introduction interpreted. The extensive field exposure of th 
National Bureau of Standards provide a considerabl« 
y of tl field corrosion data obtamed by volume of data that can be used fo his pu pos 
Ol il 4 eC: ors ras invariably results The orrosion cell ised 10] this Investigator sa 
ision that nos soil property is a good modification of the cell deseribed by Denison It 
yt ~ CoTTOSs ward iron ol ster lena tself Oo bette! pro Dility | Tl the oO} nal 
1] 1] f 
I normal corre 0 ron and steel nae cell and Will remaim opera ‘ 0 Lt lone’ T ) ot 
n electroly pl Ol ion. several accel mie Although ele ea men ement ean he 
oratory and tf based on electrical made that account fo Wwe loss of st | expose 
~ mVve 1) ‘ that puUrpo! Lo SOLS ( ht lo ( ‘ MS1OT) ey 
: , 
( Co S1\ 0 kw ny | ha 1) ius ( i} Nl ‘ 1) ‘ » Ne 
hat tl of such tests in effect rhage rr neti oO a wv the el 
( 50 ( s Whereas measure rosior te 
s STIVILY ( ong a pipe lime might 
v we th nee of corrosion on 2. Experimental Methods 
of that p ‘ ~! 30 e that 
’ th OrTosiol . e 
eneral co etwer he co 2.1. Description of the Corrosion Cell 
meta expos ) ent souls ame the 
t\ The ecorrosiol erie oO s compol ts ar 
2 des boratorv corrosion shown hie Ihe ‘ 0 sts es tialiV Of a 
h an el ' ) re was developed ched on oO eel anor \ Tit , elioh 
1 ery ) V« stee] electrodes periorates thos B I i ol thy rT) materia 
by soil as the rivte The maximum both provided vitl | connectior ind 
sures betwee! | electrodes over a 2 separated DV soil of ( Hl ‘ ost econtel hie 
od Was <4 is an dex of soil corrosivit\ soil is conta ed mn vo | e evimnders. ( Other 
ition Was ol i ( pbetweeh laboratory parts ol hie cel thn ( oO hole 1! oO rl sac 
ita exe tf } ous of relatively h of) ind thus ms ‘ ood CO " ti | ‘ ode al 
\ In such soils, | ise of local action the two stainless ee! creel 1) ame | the rubbe 
cht losses w ore significant than the stopper, H, and e rubber band, G I} con 
nts as an mae ( ) COTTOSIVILN ponents are jomned with aspl il | wl h also serves 
. ] | } 
Ewing, and Di 301 1] have shown by as ad sewiel (eration o ir ithode is accomplished 
osure tests that the corrosion rate of 1ron and by adjusting the moisture content ol tha ou, and 
sous Varies trom an almost constant rat with access of air to the anode Is rest ed by mecha cally 
some soils to virtual cessation within a few WOrkKIng the laver ol soul n contact \ th t why ‘ hy 
other soils Such information ts of great results in a difference of potential bet el the 
those concerned with underground construc electrodes. Details concerning the preparation of 
he time required to conduct such tests in the | the cell are given in the append 
owever, 18 long, and the required information 
be obtained. m many cases, when it 13 most 2.2. Use of the Corrosion Cell 
\ suitable laboratory test using the soil m _ : 
{] The behavior of the corrosion cells in refleeting 
on and operating overas ifficiently long period e 
' { the effect of time on corrosion of ron o steel is 
anv change in the rate of corrosion of iron Catteni t , : TY ' 
} rr Liustraced \ ngrure é ie cel eurrents yer 
| with time is most desirable The data from | i 
} ] measured With & Zero-resistance m immmeter wit 
Lest should be correlated with known field- ’ , t] 1 | —_ Th 
: ou vermitting 1@ cells to be on open cir ‘ 
data so that the results can be correctly . —_——— 
: currents for each soil are the average values for four 
i cells It will be noted from the table in the figure 
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that the ele trode welght losses indicate the 
7 order ol corrosivity ol the soils whereas the i 
Vv 
F vs cell currents in soils 62 and 26 are in revers: 
because ol the soil resistivities However ba 
=—_ ‘ 
Al Deteit 2 ieee field data of Denison and Romanoff 6) the ew 
{ 2 ; 
reflect the effect of time on the rate of COrro 
c Puddied Soil ferrous materials in the three soils. The fiek 
Detail |) + reveal a fairly constant rate of corrosion of wi 
Seneted Oat materials exposed to soils 64 and 62 over a pe 
F 14 years and a stifling effect in well-aerate: 
B( Detail 2) similar to soil 26 It will be noted, figure 2, t 
soil 26 the cell current gradually decreases, W 
Drill & Tap-.250° Deep in soils 64 and 62 after the initial decrea 
meme Seren ~ ee currents are relatively steady 
" ae OC) OX notes.250" Dio As a result of this and other experiments, 61 
- ‘| 2 (QO CO : was chosen as the minimum time requir 
ey? ei 3 lOO0GD)O measuring soil corrosion Although the labo 
2 | | ae 7 CE OO Oj/rosius-o:.344" results are based on the weight lost by th 
* L --* i - . +S. . . 
= 2 OAC Rodius- b=.688 electrodes the cell current serves to rey 
. ~ 1766" * 020 X quality ota cell If i divers cell indicat 
— 1-750 oo incorrect polarity or relatively low current, « 
a particular soil vroup weeks after assembly 
DETAIL | DETAIL 2 that cell should — be replaced Six months 
2"K3/16" Lucite Tubing Electrode A-125"th.t.010 assembly the corrosion cells are taken apart al 
One Side-Finished Surface electrodes cleaned by methods 7 already (Lest 
Electrode B- Some As A, The average combined loss of metal on the 
Except With Holes and cathodes for the four cells is used as the e 
uy index 
3. Comparison of the Corrosion of Steel 
, Cast Iron 
3.1. Weight Loss in the Laboratory Tests 
6 
Although the laboratory corrosion test was pla 
primarily to show the relative effects of Varo 
6 5 
Pa on ferrous materials, to be shown later, the sens 
~ of the test can be better evaluated if it is pos 
= 4} to at tect ithe rences in thre corrodibility ot com 
- \ " ang used plain ferrous materials when exposed to a 
= 3} ee 4 . - a = 2s soil environment To that end, corrosion cells 
_ ° ca > assembled with both steel and cast-iron elec 
4 | - 
Ww 
© 2} PAR. Pp 
ot} . 
ie 
8 e—e- o~ 
. o—-¢ 
——— ! eodinniadlian SiON eee ee 
.) 20 40 60 80 100 120 140 160 180 ( 
TIME, days 
“- h 
hy re 2 } j ¢ it lay adot yak ( ‘ 
. Far * " I 
H black cla Ar 1s 
M ia la al M ilwauk W * 
Re r f ( f 2 | Ontar Rochester, N. ¥ ‘ 
. il ‘ t loar Ka City, M 
‘ 4 Li abeth, NJ 
Ha wl i Lock Raven, M s 
Ohm ‘ n ‘ Lake Charles clay El V ex 4" 
a i 62 ns M t harkey clay New Orlear La ) 4 
2 122 f Intermediate #2 jUehanna clay Meridian, Mi ‘ 6, 922 
. ” Patt 7 Least 1 Docas clay Cholame, Calif 5 62 
Chino silt loar Wilmington, Calif 7.2 148 
Calculated fr the area under a through the 
Combined weight {the anode a hode ex} 180 days, in ounces es im ber al fe t im ber as specified the NBS 
er square { udic 
Order shown by t! ld-exy te a erat f G F, fair; P, poor; VI poor 


330 











i i> 
National Bureau of Standards and known to of 2 de Lavaud specimens i northern pit cast iron 
he range ot corros t\ Che location of the and 1 southern pit cast-iron specime! Data n 

tes and some ot the so properties are given in table > pertamimg, to soils 56 H2 eat 4 vere 

niong with other so to be discussed later reported D Denison and Romanoff [6 

results of the laboratory tests are given in In all soils (table 3) except soil 62, the teel spee 

2. Inall of the soils, with the exception ol soil 5 mens lost less weight than did the cast-iron spe¢ 
electrodes were more resistant to corrosion mens he statistical significance of the difference 
ere the cast-iron electrodes. The significances n weight loss. as in the case of the laboratory com 
lifference indicated in the table is based on parison s based on thr test and, except in ous 62 
indard f-test [8 t will be noted that the 13) and 64 hows substantially the same result 
corrosion resistance exhibited by cast iron For soils 43 and 62 thre aboratorv and tield com 
n soil 5 may be attributed to chance parisons of the corrosion of steel and cast tron are 
neonclusive: and in the case of soil 64. the field data 

Comparison of Laboratory and Field Data | are so variable that anv significant difference be 

tween he steel and i ron is obs« ‘ | oil 
ymparison Of laboratory and field data pet 64. the standard errors appiving to the laboratory 
: _ og Meg omnr - pti — ron “ms data (table 2) are relatively small Pherefore, the 
mea 2 SOU Sampecs are _ laboratory res ilts in soil 64 are to be taken a Heme 
ve of the soil in the trenches from which the 
specimens were removed The field data are the more significant 

n table 3. For all of the soils, except 56 The preponderance of negative values in column 
nd 64, the data are based on weight losses | 8 of table 3 is not to be considered as indicating the 
ted by Logan [7] occurring on pipe specimens relative behavior of the materials in soil enerally, 
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4. Measurement of Soil Corrosivity 


Correlation Between Laboratory and Field 
Weight Losses 


In order to further evaluate the laboratory corro- 


41 


relative corrosiveness of 
materials, the 
mean we ight losses of steel electrodes afte exposure 


sion test as a measure of the 


toward commonly used ferrous 


soils 


for 6 months to the soils shown in table l were 
correlated with the average weight losses of other 
specimens exposed to the test sites for 10 years 


aled S in number 
materials 
iron 


The At ene h test site tot 
and consisted of 2 of each of the 
st¢ el 


Sper Hens 
follow Wie 
hand-puddled wrought me- 


and open hearth 


Low carbon 
chanically puddled wrought iron 
ron 

In developmg the laboratory Lest thre National 
Bureau of Standards having the 
cooperation of the Waterways Experiment Station, 
Army That 
followed the same procedure in measuring the corro- 


The soil 


Experiment 


was fortunate In 


Corps of Engineers, | » laboratory 


siveness of soils as described in section 2 
tested by the Waterwavs 
Station included 9 of the 10 studied by the 
National Bureau of Standards and in addition 5 
soils removed from other National Bureau of Stand- 


sample Ss 


soils 


ards test sites The 9 soils tested by both labora 
tories were not taken from the same samples but 
were removed from different parts of the trench at 
the 9 sites Both laboratories used steel cell ele 


trodes of the same chemical composition and dimen 


sions The composition and dimensions are given in 
table The data obtained by both laboratories 
and the National Bureau of Standards field data are 
tabulated in table t Correlation between the 
laboratory and field weight losses is based on the 


The coefficient of 
National Bureau of Stand- 
ards laboratory and field data is 0.84 The coefhi- 
client between the Waterwavs Experiment Station 


text of Croxton and Cowds nis 
correlation between the 


PARLE t ( 
\ or Bure Star 
‘ ‘ 
W 
‘ 
} 
. ; 
‘ i 
\/ Vl 
ul ‘ 
ae 
" « ’ 
‘ 4 
‘ % 
wu 
4s 
; M4 
~ % 
‘ ‘ ‘ 
. 
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data and the National Bureau of Standards 
data is 0.6, but by eliminating soil 64 the coeff 
becomes 0 S7 

Because the souls listed in table 4 cover the ral 
corrosivity of soils peculiar to the United States 
data of both laboratories are used to set up the ¢ 
tion of best fit [9], presumably a linear relatior 
between the laboratory and field weight losses, w 
it Is again presumed will be applicable to any 
The Stam 


losses were at 


whatever its degree of corrosivity 
errors associated with the weight 
the same for both laboratories ex« ept for the W 
ways Experiment Station data pertaining to so 
in which the standard error was considerably la 
than for all other data As the 
64 adversely affected the correlation 
Waterways Experiment Station data and thi 
data, the Waterways Experiment Station data 
Thus by uvera 


inclusion of 
betwee! 


taining to this soil are omitted 


the weight losses in each of 8 soils (4, 5, 25, 26 


56, 62. table 4 tested by the two laboratories 
using the individual laboratory weight losses in 
7 other soils s. Ld. 1] 95. 61 64. 65 table 4 


correlation between the combined laboratory we 
losses and the National Bureau of Standards we 
losses in the field results in a correlation coefficie 
O.S5 All correlation coefficients wert 
Based on the 


re lating the 


correcte« 
size of sample correlation coeffi 
of O83. the 


field data is 


equation laboratory 


i S.45u l 





the anticipated weight loss on 


where W is 
commonly used ferrous specimens after 10 year 
field exposure 


wis the mean combined weight loss in 6 months 


expressed Im ounces per square 


the anode and cathode of corrosion cells, expresses 
The standard error of « 
mate associated with this relationship ts t.1 07 

The the Natio 


Bureau of Standards specimens based on eq (1 


ounces per square foot 


calculated weight losses on 


tabulated in table 4, along with the actual we 
losses Relatively large discrepancies appeal 
tween the actual and calculated values im soils 
f1, and 61 In the case of soils 41 and 61 


discrepancies might be attributed to the laboraté 
test Soil 5, however, is fine in texture, very de! 
and poorly aerated, and also possesses other char 
teristics which Denison and Ewing 
with very corrosive soils. Also, a 
early information published by Logan, Ewing 

Yeomans |{11] revealed that iron pipe removed f1 
soil 5 was severely corroded Probably the condit ( 
within the trench from which the field 
were removed were not typical of the general 
structure Therefore, for soil 5 at least, the labe 
tory data considered a better index of 
potential soil corrosivity than are the field data 


10] associat 
reference to si 


specim S 


4.2. Correlation Between Weight Loss and 
Maximum Penetration 


The value of a laboratory corrosion test wouk 
immeasurably increased if it could be 
predict the maximum depth of pitting under no 


used 














Laboratory weight losses might bye 


exposure 
terms of pit depths if it could be shown 


rpreted in 
correlation loss and 
after 


2 specimens 


existed between weight 
on terrous 


this possibility, 


mum penetration specimens 


exposure To test 
$ commonly 
previously described, exposed for 10 
s National Bureau of Standards 
Following the stat stical procedure already 
the of the deepest pits 


Thes¢ Spec lens were 


ach of ised wrought ferrous mate- 


Vears at 
test sites were 
sen 
depths 

correlated with the 


read to, average 


average 


ht losses The coefheient of correlation was 
ited to hye 0.74% he eq imntion based on this 
ation, relating weight losses and pit depths for 


alter 10 vears ol 


8 soils xposure 


P b.2 6 ” 


Py is the maximum pit depth in mils on a 
i O41 and WW s the weight loss expressed 
nees per square Toot The standard error of 
ite pertain ng to P? n this equation ts IS 
Because of the mber and variety of soils 
ed in the derivation of eq (2), the relationship 
be considered as beimg generally, applicable 
» wrought ferrous materials exposed to any soil 
\ tional evidences of some correlation between 
ht loss and pitting on wrought and also cast 
Is Spec mens s ate ated Dy the Studies of 

l) son and Romanoff |6 
By combining eq (1) an >), the maximum antic! 
it depth at LO vea s related to the labora 

it loss for 6 months by the relation 


Oh an CN 
lor 
maximum 


eq | 
actual pit 


field specimens and also the 


as expressed 


ited values based on eq (3 The lack of rea 

I ood agreeme! n soil 5 1s to be expected 
se of the facts pre ous pointed out concern 
his soul 


5. Adjustment for Time and Exposed Area 


It has been shown that b means of the laboratory 


{ ble estimate can be made of 


rosion test a resona 


veight loss and the ave race cle pth of the de pest 


that 


posed m the 


and steel specimens ex 
the prac 


will occur on tron 
field for 10 vei 
value of the test 
ssible to predict similar information as applied to 
| exposure other 


Obviously 
enhanced if it 


Ls 


1 


would be were 


Cl exposed areas and periods ot 


than 10 vears 
Martin [12] has shown that 
th veight loss of iron and ste« 


may be expressed by the empirical 


the relation between 
| in soils and the dura- 
on Of exposure 
ation 


Wang 7” 


is the 


and wu is 


time 7 k’ 


tire 


where Wis the weight loss at any 
weight loss on a unit area for a unit 
a constant peculiar to the soil 
Similarly, Logan, Ewing Denison |4] 
shown that the average depth of the deepest pits on 
be expressed by the empirical 


ana have 


iron and steel may 
equation 


P=—fI"A 


where P is the average depth of the deepest pits at 


any time 7 on the exposed area A, and & is the 
average value of the depth of the deepest pits on 
specimens of unit area for a unit time. The ex 
ponents 7 and @ are constants depending on the 
aeration and drainage of the soil, and @ beime rela 


tively fixed in all soils 
Kor a practi al application of these equations to 


soil values must be selected for the 


any auppropriat 
exponents Logan, Ewing, and Denison tabulated 
values for exponent a based on exposure tests of 


characterized 


wrought ferrous specimens In 17 soils con 


ided that the value 


otaco ild not 


by any particular soil property As the values re 
ported did not vary over too wide a range, the writer 
believes that SOT averavt value would bye appro 
priate for all soils. By eliminating all soils in which 


tabulated ior d eX 


for the 


ated standard errors 


thre 


eeeded remaining 


20 percent the values of a 


s1 souls ranged between 0.082 and 0.241. averaging 
0.15, with a standard error of 0.019 

Logan Kwing and Dension divided the 47 soils 
into 4 Yroups based on aeration and dramage and 
calculated the average value of the exponent for 
each rroup These values ol and the associated 
standard errors are shown in table 5 Average val 
ies for uw and standard errors based on 14-vear ex 


rials b Denison 


posure tests of wrought ferrous mate 
and Romanoff [6] are included in the same tabk 
It will be observed that the values of both ind wu 
bear an inverse relation to the degree of soil aeration 
| RLI 1 i , 
° ° , 
h Hi. Log Pp jf 4 
ASTM j 
. MA 
| } “4.4 I 





Also. if the values for ana ' ompared, giving 
consideration to the standard errors, the mean values 
ior 7 and uOovVve rlap CAC pt lor yroup II] Beca ist 


the values for uw are based on relatively few soils as 


compared with the values shown for the later are 


taken to be the more significant For this reason 


has already been shown concerning 
weight 
believes that the values for 7 
applied to 7’ in both the 


\ccordingly 


of what 
the correlation 
penetration the 
ith table ) 
weight loss and pit-depth equations 


be« nmus¢ 


between loss and maximum 


write! 


shown cuh bye 


this is proposed in table 6, together with the dis 
tingushing soil properties pertaining to the four 
froups 

Having obtained from the laborator test a value 
for the average depth of the deepest pit at LO vears 
Pi on 0.4 ft? (eq 3), the value of pitting for a time 
T’) other than 10 vears and for an area (4 reatel 
than 0.4 ft? is calculated from the equatior 


ij A 


Pp Po( 39) (oa) { 


It has 


Ly hnison 


been pomted out by Lo an Kwing 
that substituted for 7 and A 
chesen indiscriminately This is particu 
substituting for A, because P 
\ reasonabl extrapolation based 
table 4 might be 


ana 
the values 
cannot be 
larly true when 
applies to 0.4 Tt 
considered 


”~0 Itt of 


and 56 
Assuming an area va ot 20 Tt 


on soils 55 


“upproxrx 


’ in steel pipe and an Indergrouna-ex posure 
period (7°) of 40 years, the maximum predicted pit 
depths calculated from eq (4) would be 112 mils, and 
660 mils for soils 55 and 56, respective In thes 
calculations, values of 7 0.19 and 0.68 were 
chosen for soils 55 and 56, respectivel the values 


being selected from table 6 on the basis of the soil 
Thus n the case of 


In pit ce pth after LO 


properties found in table | 
soil 55 the increas ars is due 
chiefly 
resulting from soil 56 are 
effect of time 
Weight losses at 


can be 


to the increased area, whereas the dec per pits 


to the 


attributable mainly 
based on the laboratory 
test (eq | similarly and 
weight loss (Vl pred ted for periods T’) other than 
The weight loss becomes 


LO vears 
aayjpusted values of 


10 vears 


/ 


(iG)? 5 


where wu is selected from table 6 


The effect of 


uM W 


on the soil 


based 


loss is 
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to be 


This equation is parte 
the 


average 


usSsunLed neat 


useful in) calculating 


Thus the 


average penetration 


conversion penetration mn 


bee oOmMmes 


P=1.53 W 


where H 
eq 5 The 


s the weight loss in ounces per square 
predicted average 
steel surface exposed to soils 55 and 56 
> mils and 81 


isel il whe re 


penetration or 
table } 
mils, respect 
} f 

i 


Oss oO 


10 vears would be 


Such information is strenet 


thre primary consideration 


6. Summary 


\ mocdthication of the ey nison soil COTTOSLON (¢ 
The 


both steel and cast ron electrodes expose: 


described cells were assembled in quad 
with 
varving degrees of corrosivity Be 
ferrous 
function of time all measurements are 
weight loss of the 


eXpos ire 


10 soils of 


the eorrosion. rate ol metals in soils 


base a on 
months 


electrodes niter 1) 





In 9 of the soils. the steel elec trodes lost less we 
than the cast-iron electrodes Statistically 
pressed the difference in weight loss in 5 of tl 
soils being aue to chance was less than 1 per 
and in other soils less than 4 percent 
laboratory results compare favorably with 12 
field exposure tests 

Good correlation was also obtained betwee 


the corrosion cells and the we 


weight losses of 
losses occurring in the field on other wrought fer 
specimens exposed for 10 vears to the 
Based on 10-vear field exposure of wrought mat: 

at 5S test sites throughout the | nited States, a { 
found to 


penetration 


same so 


correlation exist between weight 
and maximum 
latter correlation, it 


used to pred t 


was 
$y making use of 
is shown how laboratory wi 
probable maxim 


Ipon field eXpos 


losses can be 
pit depths on ferrous materials 
alter 10 vears 


$y soutilizing relating weight loss 
ferrous field specimens to time of exposure, and 
average depth of the maximum pits on the specin 
to time and to exposed area, values of weight loss 
the end of 10 vears predicted by 


laboratory test may be extrapolated to any des 


equations 


pitting at 


area and period of exposure within reasonable lin 
the 
adjustment of the data are given. 


Suggested values of exponents required 
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8. Appendix 


Preparation of the Corrosion Cell 
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The 
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Effect of Temperature on the Electrical Resistance and 
Voltage Departures (Errors) of Glass Electrodes, and 
Upon the Hygroscopicity of Glass 

Donald Hubbard 


fe perature uy e ¢ 
hed prepared I! et 
" i { ng O15; glasse { 
‘ ‘ he s how na ele I i 
int wwe aeparture tron ‘ etica 
4 4 4 fa indicators of lrog 
e ¢ ie epared from Pyre ass ga 
sta i i vy as 12 megohms The re 
vgroscopicity gave a pseudo-correspondenc 
rode ( ng O15 gave voltage depar 
a i e ot the electrode 
1. Introduction 
(‘onsiderable confusion exists in the general con- 


as to what part the electrical resistance of a 
ss electrode plays in its voltage response to hydro 
of aqueous the voltage 
errors) that appear high alkalinities 
n the “superacid”’ region, and particularly in the 
of from many to follow 
simple dictates of the Nernst Ak 
00198 TApH For example the difference in 
rformance shown between electrodes from Pyrex 
fig. | 
resumed to be related to a difference in 
of the glass membranes forming the effec- 
portion ol the electrode With the electrometer 
bes available Ml 10 
amp trode having a resistance as 
as 1,009 megohms would exhibit an IR drop 
ror) ranging from 107° to 107° v {1 17 and 
Such most 
cal pH measurements can be disregarded It 
that exhibiting large 
ge departures, pH errors, but having resistances 
than 1,000 full 
electrical 


on activity solutions, 


mart at 


ires 


electrodes glasses 


' 
ire 


equation, 


and a satisfactory pH-responsive glass 


> 41) 
SOTLE how 


sistance 


requiring gt irrents of only 
() a glass ele 


pages 


, to 92 errors are small, and for 


MVIOUs class electrodes 


megohms 
not 


are 


failing to show 
ise of their high 


oH esponse beca 


stance nor because of a lack of sensitivity of the 
able indicating electrometers 

he present investigation was undertaken to as 

tain what individual parts the hvgroscopicity 

sorption) of the glass and the electrical resist 

of the glass membrane play in the over-all 

re of the ability, or lack of ability, of electrodes 


t as indicators of hvdrogen 


Various classes to wa 


tivit\ ol aqueous sol itions ove! the temperat ine 
20° C to near the boiling point of water 
2. Experimental Procedures 
| hyvgroscopicit Values were obtained in the 
aratus deseribed in previous public ations [2 


ical re tance and pH t ( is Det 
pe glass¢ Pyrex 7740 ercia f 
inte liate and high hve s pieit f 
re pa | ( ( le al oO 
and ere Ta ‘ Toma BS 
aC ia i> For « Arthy 
gnificant pH respor en whe 
of ele les from glasse rit ediate 
f pif ( ‘ \ resista { wi i 
Ir na | iV be at ed 1 


The elevated temperatures 
carried out by immersing the humidity chamber in a 
thermostated water bath Although no mechanical 
difficulties encountered, this for 


obtaining hygroscopicity was not entirely satisfac 


experiments al were 


were procedure 


torv, because of condensation on the walls 
of the chamber at the elevated temperatures for the 
high humidity (approximately 98 percent 
for the experiments. Howeve1 
to be qualitatively acceptable 


those 


CXCECSSIV¢ 


required 
the results appeared 
The hvgroscopi it\ 
values obtained for 1l-hour ex 


reported nr 


pos res 
The 


electrode 


of experimental 
the voltage change that 
upp vared when the cell consisting of the « xperimental 
electrode and a reference electrode was transferred 
from a buffer solution at pH 4.1 to a buffer at pH 8.2 
The reference electrode emploved was « ither a hydro 
at the 
experimental 


voltage performance any 


was determined by 


ren electrode or a 
elevated 


glass electrode calibrated 


Thus if the 


temperatures 


clectrode gave full pli response no change in voltage 
appeared in the two buffers at any temperatures 
Typical examples of the data obtained by this pro 


COMMcre ial 
7740 


illustrated | for a 
Pyre x 
pli interval was chosen because it 


of thre 


eliminate d or 


cedure are in tigure 


electrode class 
The 
region of 


I idied 


ana om 


pre pared from 
above 


uniform chemical durability classes 
Thus, voltage shifts are 


minimized, which m ohit be credited to a change in 
chemical durability [2 The calibrated mercury 
filled glass electrode [4] was preferred over the 


electrode becau 


ol 


hvdrogen elec trode as a reference 


operational simplicity. and ¢ 


0 usé temperat iné 
control by eliminating the hy lrogen-yvas cham 
Furthermore, with both the reference and expert 
mental electrodes filled ith mereury any small 
thermoelectri voltages venerated at the Pt Hy 
junction would be compensated by the svmmetry of 
the svstem All electromotive-force measurements 


were th a Beekman pil meter Mmborator 


Triad W 
model (; 
The 


elec trodes 


th 


the 


experime ntal gl 
different 


for 
rT 


resistance Values iss 


were followed tempera 
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tures with an RCA Senior Voltohmyst TW) 


which imposes 1.5 v across the glass mem 
while the bulb of the electrode was immers 
Britton universal buffer at pH 4.1. The ek 


connections for these measurements were estab 
OV means of platinum leads, one of whicl 
dipped into the buffer solution and the ot] 
the mercury filling the glass electrode bull 
number of parallel observations were mad 
General Radio type 544 BSS megohm bridg 
ating with 100-v drop across the terminals, to 
sure that neithe polarization nor the voltage 
erated by the battery consisting of the two 
glass electrode, and platinum wire was i 
spurious resistance values to be reported [5 
higher voltage of the megohm bridge might « 
polarization resistance at the glass surface o1 
puncture the glass membrane, whereas the 
operation voltage of the Voltohmyst might 
seriously affected by counter voltage veneratec 
the battery consisting of the glass electrode an 
external platinum-wire connection Because thi 
instruments indicated approximately the sam 
sistance, it seemed reasonably sure that the eff 
of polarization and interfering voltages wer 


lara 
3. Results and Discussions 


3.1. Electrodes Prepared from Pyrex Glass 7740; A 
Glass of Low Hygroscopicity 


Figure 2 and table 1 show the effect of temperat 
upon the electrical resistance and pH response 
electrodes prepared from the sample of Pvrex 774 


chemical-wart glass and the hvgroscopicity ol 





glass over the same temperature range Phe 
sistance curves for these three electrodes ar 
accord with the long known fact that the elect 
resistivity of glass decreases rapidly as the temper 
ture is increased and follows in a general wa) 


( xponential equation 


R= Be* 
where 
R=the resistivity 
B—=a constant, 
A=a constant involving the activation energ 


r temperature, 


é€==2.71828, the limit of the series (1+-1/7 
n=ea 6] 
At times the resistance curve obtained experiment 


may depart conspicuously from this shape; typ 


examples obtained on an electrode are illustrated 
A and B in figure 3. Such S curves can be rea 
accounted for qualitatively by electrical surf 
under conditions in which the resistance 


great 


leakage 
the electrode bulb becomes increasingly 
comparison to the surface resistance of the 
tube The total resistance, R, of the two 
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elec trode b ilb and glass surtace in paralle! 
nsert whose separate resistances are ana 


Ss oiven by the relation 


l R l | Ol R 


of the glass stem with appropriate Water re pellants 


Sule h as paraffin organ chloros lame Ss, or cotninercila 


\ becomes so very large as to approximate the Desicote |6 
sistance of an open circuit R approaches ry as a The voltage de parture characteristics of electrode 
ting value | as affected by temperature are also plotted im figure 
2 Krom this curve it is evident that there was no 
1/R=1] | or R improvement in the pH response of this Pyrex elec- 
trode although the resistance fell to such low values 
i hese curves are readilv brought back into con as 99,35, 19, and 12 megohms In other words, the 
mity by cleaning, drying, or treating the surface | low resistance was not accompanied by a corresponds 
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ing improvement in pH respons Electrodes 2 and 
> were checked at st) nna YO ( hr were also found 
to be unresponsive to the pH of the buffer solutions 


the 
the 


ve lop 1m 


electrical resistances were 


in pH response of this rlass 
Pyre Xx elec trodes would hye 


Obviously | 
determinant three 


expected to «ae 


proved response at different temperatures mm accord 
with their res Stance curves Llowe ver, thev are all 
useless as indicators of hvdrogen-ion activity Only 


upon a study ol the hvgroscopi tv curve obtained 
on Pyrex glass does one get a plausible reason for 
thy complete madifference ol these electrodes to the 


lowered resistance at the elevated temperatures (hg 


2 The water sorbed by the glass has not changed 
appreciably with the increas nh temperature and 
does not at anv te mperature approach a value com 
mensurate with an effective electrode glass 5. 7 


Certainly thre lowered resistance can mot De attributed 
fo water penetrating the glass As a result the 
voltage cle parture of the electrode mv pH ay 


This inh 


comes greater with increased temperature 

creased departure is largely the result of the temper 
ature effects inherent in the derivation of the Nernst 
equation Based on previous evidence the elec 


trodes pre pared from this non VoeTos op cl iss prob- 
regard- 


blow n 


ably would have no satisfactory pli response 


less ol how thin thre cr| membranes were 


iss 


oo. 


3.2. Electrodes Prepared from a Glass of 
Intermediate Hygroscopicity 

Figure t effect of 

resistance and pH characteristics 


the 


Iwo glass elec - 


shows the temperature on 


f 
ol 


trodes prepared from a commer ial soft lass ot mte! 





mediate hvgroscopicity. The performance ol these 
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Figure 4. Effect of temperature upon t al resistar 
hygroscop and pH response « ectrodes prepared fron 
a glass of intermediate hygrosco; 


major 


crore ral Cones 


nearly follows the 
of the part that resistance should play in detern 


electrodes more 


the suitability of a glass fer pH measuren 
Electrode 1 at 20° C was verv sluggish in d; 
toward a steady voltage whereas at 30 and 4 
the performance was still sluggish, but mas 
improved Above 50° C the resistance of the 


fell rapidly below 1,000 megohms, and the eleet 
performed in an acceptable manner for high-ten 
ature pH measurements. As the chemical durab 
of this to that ol Corning { 
the life of electrodes prepared fror 
when used for high-temperature pH measurem: 
V create! 1, 10, 11 Ke 
trode 2 was blown with a bulb considerably thi 
than elec trode | \t >) C" the voltage depart 
was large and verv unsteady This electrode sho 

inserviceab 


rlass is far superior 


expectancy 





would bye corre sponding 


marked improvement, but was still 
10) 50° © The voltage depart ires 
at 70° C, but the 
still sluggish. Beyond 80° C 
bulb had fallen well below 
over-all performance of the electrode was acce pts 
As the pH response of thes elec trodes mproved 
of the pH 


resistance values would apparel 


and 
oH) 


were sl! 
ele etrode 
the 


1.000 


and perio mance 
resistance ol 


vy yams 


Thi¢ ane 


decrease im resistance, a correlation 


sponse W ith the 


be acceptable, but not entirely adequate, as will 
seen from the discussion and evaluation of the « 
in figure 7 However, the hvgrose« 
citv-temperature curve is entirely different from t] 
ol Pyrex and probably s equally 
accounting for the improved 


trode Compare he 


shape oft thre 


satislactory 


performance ol 
hveroscopicity « ] 


glass ( le ( 
in fig. 2 and 4 
3.3. Electrodes from Corning 015 


Figure 5 shows the electrical resistance and volta 
performance of three electrodes prepared fron 
single sample of Corning O15 glass Electrods 
had low resistance values at all temperatures, rang 
from 115 megohms at 20° C to 0.59 megohm at 90° ¢ 
Having prepared from a of relat 
high hy groscopicity and of uniform chemical du 
bility over the pH imterval chosen, this thin-wa 
satisfactory performance (appro 
mating the theoretical) over the temperature rar 
from 20 90° C Electrode 2 was blown from tl! 
same length of glass tubing 
much heavier wall thickness 
for this elec trode fell below 
and showed the usual rapid decline with increas 
temperature. However, this thicker electrode « 
hibited appreciable voltage departure at all tempe! 
tures up 70° C, beyond which temperature 
paralleled the performance of the reference elec trod 
Electrode 3, with a wall thickness in excess of 1 m 
still had values lower than many 
fac tory electrodes, but fave voltage departures oO 
the entire temperature range from 20° to 90° | 
This performance is in accord with the 
knowledge that every glass has a critical thickness 


glass 


hee Nn 


electrode cruve 
to 
as electrode 1, 

The 


500 megohms at 20 


but W 


resistance ¢ 


to 


sal 


resistance 


gen 


which voltage departures appear [8]. Such failu 
of an electrode to respond to hvdrogen-ion acti 


certainly cannot be attributed to electrical resistai 
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electrode because at all temperatures the 
stance is at values far below that at which the IR 
»across the membrane wo ild be as creat as | m\ 
Possib!y proton o1 H.O” resistance due to the thick 
ss Of} the membrane leave k ot watel penetration ol 
he glass, or a combimation ol both vive a more 
eptable concept It seems probable that the 


adrovwen-1on-actrvity 


h equilibi um does 
become throughout the glass even 


ho hh the nherent hygroscopicity of the glass is 

( th (Tig ) 
In order to further emphasize the lack of universal 
I espondencs between voltage departures of glass 


rodes and electrical resistance of the glass mem 
determimed 
extended pH rave 


temperatures trom 20 


esistance values were 
of Corning 0 
rie H.SO, to pH 2) lor 
MU) ( the 
mventional voltage ce 


The I lor ay 


rods 5 over an 


oOo 


results considered in relation to au 


and 


vmarture curve obtaimmed a 


temperature ovel the same pH range (table 2 


he 


? 
‘ 







































that the 


and superacid | 


These results show 
the high alkaline 
iot attributable to excessive ly high electrical resist 


voltage depal 


ires im unges are 


they associated with marked changes 


The 


hor are 


n electrical resistance causes have beeu dis 


issed elsewhere n connection with the wate 
tivity of the solutions chemical durability of 
he class : heterogeneous equilibria ut the class 
lectrods sol ition nteriace 7 and eq il brative 
esponse to ons other than hvdrogen |], 16 
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| 
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onl _ Pyrex -» } 
- | ay 
a ) | 
> « / / 
e Z I ( 
. | 4 . 
® 
E wo} | 
x j 
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4 | ( 
w | i 
= } 4 j 
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3.4. Hygroscopicity versus Temperature 2. Membrane walls greater than a limiting 
ness value, which differs for every glass [S 
Although the hvygroscopicity curves for the three » Poot chemical durability of th class 2 
glasses reported must bye classified “as bemg ot un- 1. Hich concentration of electrolvtes 133 14 


certain reliability, they do warrant a brief evaluation 
If adsorption were the major fa: tributing to 
the water sorbed, the hygroscopicity curves should 
with temperature [12 As th 
curves for two of the higures 4 show a 


with temperature, it is obvious that 


tor col 
decrease Increased 
classes and 5 
liberal increase 
predommating tactors 


absorption and reaction are 


4. Summary 


Figure 7, in which the voltage ires (errors 
are plotted as the dependent variable and resistance 
as the determinant, illustrates the 
played by electrical resistance in the over-all picture 
of the pH response of glass ele Not only do 
the pH response-resistance curves for electrodes from 
different coincide for 
trodes from the same length of tubing (Corning 015 
and soft-glass tubing) fail to exhibit similar character- 
istics. From the present investigation the part 
played by electrical determining the 
effectiveness of a glass for as glass electrodes 
seems to be largely superficial. For exampk 
trodes prepared from glasses of low hygroscopicity, 
such as Pyrex pH 
function even at Fur- 
thermore, electrodes from glasses ol high hy Frosco- 
picity and low resistance, such as Corning 015, fail to 
give full response to hydrogen-ion activity if the 
membrane thickness is excessive From the data 
on the above two plus the results obtained 
on a glass of intermediate hygroscopicity, it appears 
that 1,000 are 
desirable, but certainly not sufficient to bring about 
staisfactory pH response. 

Voltage departures from the theoretical pH _ re- 
sponse can appear for glass electrodes of low 
trical resistance for the following reasons 

1. Low hygroscopicity of the glass from which the 


ll | 


electrode is prepared [2, 7 


cle part 
incongruous part 
trodes 


glasses not but curves elec- 


resistance m 
LIs¢ 
elec 
a staisfactory 


not develop 


very 


do 


low resistance values 


classes 


resistance \ alues below megohms 


ele - 
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5. Inhibiting films [9 


| 
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On the Numerical Solution of 
Parabolic Partial Differential Equations 


Gertrude Blanch 


| 
i i i 4 , 1 i 
i i I iria 
i i i 
i i 
l i : i 
i i 
i i i i WW 7 
fT i r l l 
i i! | i 
i \ i I 
i i i 
1. Definitions; Basic formulas central differen we define, 
he ( Ist ( ola nique s utiol ton ¢ ~ ‘ . 2p 
} x n nig olu n L lif x | Ju ; » 
il equation has been established, and an i Pru ‘ 
roximating function has been found that con- 
s to the solution under s table conditions ther The pomt in the s/f plane with coordinates ml 
ns the ye — providing ca — n - ab will be denoted br (aa). 
reat ent ( { ( aDDroxXxil ; ul Study rr 
aumel , ; tp] ma On , , Let uw be regular at (7 Then for & sufficiently 
cerns itself with one phase of this problem, for small there is a Taylor series in ¢ around (n 
se when the approximating Tunction Is ex 
ssed bv a differencs equation We shall furthe ] ) 
‘ : ~~ C i 
the dis ussion to a specihic type of differential ‘ a p! Of , 
ation, namely 
If.in (3) terms involving p>2 are dropped, and if 
Ou Oru L { tay () « <7 0 <i} Ou Of Is re placed by the right-hand side of (1), with 
Of Ol e O°“ Or” approximated by central differences, we get 
the well-known approximation to u 
h given initial and bo indary conditions Let us 
rodauce a lattice coverme the region, at intervals \ kf 
the x-direction and k& in the f¢-direction Let 
/ be the mesh rati It is a trivial restriction i =A \ r | / r i 
assume that a sh, where s is an integet 
NOTATION For the sake of brevity, we shall writé Similarly, if terms through p=2 are retained, we 
l if (x,t) is a lattice point obtain 
\ a 
0” 4 0’ de » 6° ae ky 
> 4 | ( 2 i2) | 
Or Ol ks Or Ol’ 
5 } 
. 4 + y 
_— (1 5) 3X?) (3 2? } 


ere g is any function under consideration 

usual convention for even 
\ . 
ke 


ermore, following the 
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where where two « xpressions for I’, are given beloy 


WJ 4 od ] a] tay °5 “+6, OU 
ay, 


Formulas (4) and (5) need modtiieation at the bound- 02 Oy \ O*u 
nu! to satisiyv given initial and bOUndary conditions i . fe” +{ ) +) 
a) of at 12 a 


to a required accuracy This modification can, in 
reneral tee Tehancie Tr the present study we shall 
assume that all required derivatives eXIst and are Where 
eontinuous Kos parabol quation his is not a WO A i ; 
( ou estriction, for let is assume that we have 0% iz U0 
enerated values ol ior a given Krom 1) and 

it is clear that these serve as boundarv values in a In the above @, and @ depend on f and its de 
subdomain for generating the set of values for th f f=0, then 6,—6,—0 Further, 0//du and d 
next ¢/ in the lattice In ln posing the continuity imply evaluation of these functions at " 
restrictions, we therefore merelv imply that regions sponding to a value of &@ intermediate betwe 
close to the given boundary vhich ma have col and Again, (6 7), (8), and (9) need moc 
ners or other discontinuities, wil e treated sepa tion in the immediate vicinit of the boundary 
rately. This in fact must usually be done in practic DEFINITION / will be said to be of ord 
either by choosing a lattice that much finer than it contains /#° as a tactor, but not / 


that required over the majo! portion of the region or 
by special approximation that are auppropriat lol 
the particular problem. Our main concern here is | 2, Stability Considerations; Bounds for thé 


l 


with the choice of the mesh ratio, A, for the majo! Errors in the Approximations 


portion of the domain where the function ts presumed 
to be regular The manner in which an error (01 Deweneen Th aleatins __ 2 : 
mN e s l ] ’ uhh upp 
atio t I i conditions is propagat : 
varia ; mn) in 7 ~~~ ir ! pre pase d mating differenc equation will be defined as 
ove t rest { : i i must, ot cours ( - } 

: ois set - domall = _ . pe ON if it is bounded for all finite ¢. independently of 
amined; but this problem is no different at th Consider th special case when fir tu () 
boundary of the given domain than at any of the ’ , , ( 

: “(zr,0 A(x) 1s defined and bounded for 
subdomains (that is, at the successive values of ¢ In 


' T) » - It will be convenient to refer to the differential eq 
the attics IS problem Will no em considered } hasic 


} unde! these spechil conditions as the fi 
ere 


Jdeneous equation Let this basi homovgenes 
Truncat ery Let It can cana } ' . 
. equation be approximated by a difference equat Ol 
iy verified that corresponding to } » ‘ 
order 2p, so that we can write 
0 
\ k] 6 SY Ss 
Ju homed ya 
where 


where the coefficients 6, and a,, are suitable consta 


\ l \ Ot | ; 
Tis } ( ) A | A) 7a It is easv to show that La i; In general 
, 


2 ie Ju j coefficients a, will be functions of X Let q be 
, uppel bound of Ala If it 1s possible to choos 
Another form for 7 is given below so that all the coefficients a, are nonnegative 


shall have, from (11 





1 OV 1 Of 
Is « } | . Oth 7b 


Corre sponding to >). and by induction on n we can establish that 
for all m,n; henes rt) is stable The condit 
, . that all a, be nonnegative its therefore a suffici 
0 = at AO ( = : ) ) though not a necessary condition * for the stab 
ns Le of rt in the case of the basic. homogene 


ou it the Inst N Analy 
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We shall refer to the range of \ for which 
asic, homogeneous equation is stable as the 
ssibl range ot A In what follows \ will 
s be chosen to vithin the admissible range 
av not guarantee that the solution f) will 
i lor arbitral OUT il conditions and tune 
f Howe I ehoice ot AW th I} th s 
isually simplifies the error analvsis for any 
boundar conditiol ven when (1 s not 
ous 
ning now > Ul pproximation formula 
by t tis clear that the admissible range of 
) \ Phe co esponding equations tol! 
n (6 Li pper bound of 7 
be an upper bound ot Of ou this uppet 
to be independent of / If o 0, eq (6 
oO 1 Wide ! f initial and boundary 
ns Since \ 1s in the admissibl 
' | 1} 
] 
| Ve: hence a ne 
) 
ma ons \ mpose modiiication of 
If cond : such that no nevatiy 
his added as a tac wr to the right-hand sick 
nequality 2). ther an be shown from thi 
/ tha i DD bound of can be tound 
| is a facto It follows thet on 2 
1s | ) 
. how co second proxima 
rmula, defined by (4 The coefficients of 
{) 2 rn ‘ ill Db POSITIVE 14) i 
roximatlo » the basic homogeneous 
mn will therefo stable for this range of \ 
\ he same inalvs is betore we can show 
approach s/ approaches ero, tol 
inge ot bo rv conditions 


the Choice of a Suitable 
Mesh Ratio 


Criteria for 


of A.A, and of the bound 
th the 


upper bound ol 


rrol isa lunctiolr 


onaditions associated W 


differential equa 
Griven an that can be 
ted in the the problem involves choos- 
and X (the within the admissibl 
is to meet requirements with the least amount of 
We shall assume that for a given approxima- 
x heme, the work ts proportional to the number 
tice points at which wu must be 
s is not strictly tru For let us define a profil 
set of values of for a fixed n 
for example, successive profiles are generated from 
eding ones on an IBM machine such as the card 
grammed then merging operations 
iv be required at the end of a profile which may 
ynsume Thus a grid of 10 points in the 


eTrTo! 
solution 


latter range 


evaluated 


j 


and all ii 


calculato 


some time 
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r-direction and LOO points in the f-direction may take 


more time to generate than a ad of 20 points in the 


r-direction and 50 points in the f-direction Neve! 


theless, the assumption that the wor Ss proportional 
to the number of lattice points is close enough to 
reality to be useful OT cours 1¢ COMPLEXIL of the 
programming must be considered; thus an appront 


mation formula of order four may take more machine 
evcles hence more tim than one of ord two 
Hlowe Vel lo! thie sume approXimMmation cheme hae 
choice of X} does not char th oul 1 wor 


he range of The number of " point thre 
region is Ns—(Nt,/hk)—Nt,/Me. As Nt, is fixed 
he wor req read tor a napproximatiol heme 
s therefore nverse! propo ona oO Al 

If the exact sol ition tor were nownh, it would 
be theoretically possible to stud the magnitude of 
the error for various chon of X and / correspondu 
to a vel pproximation schemy The pi 
solution snot easy to find Llowever, from P 4 

scleal hat ay Ipper bound which can b approx 
mated has / as a tactor We shall, therefor 
ium to choose XX and / n such a manner a to make 
/ Sinai \loreovet for both approXimation 
formulas (4) and (5 the successive terms of h?7 
nvolve / ) Lo and / »’* ?t/d We shall 
require that for all choices of A, th nterval fh b 
s fie I md | i sO) - Sil siy 
QiIMOSL ¢ vh | ) th how ons con 
rhe ne «ce il \ h differences Ss cleu hat 

mpl s ha - ( SSiveé ms of / 1) 

numerica smaller than precedu one I} 
phrase almost everywhere for tl conditiol | 
needs explanation. It may happen that in a region 
where changes sign, a lew entries ol may 
be numerically matller than corresponding entries 
in the higher difference Such a case ma ilso arise 
near critical points ly parti ir th Orne 0 
shall be satisfied b hie ! ) na 
) In practice one olte q ( that eb Ol 

For one power Len ) l Te ( ( om 
puted \ ilues is obtained from th patter ts ces 
sive differences of the entri He even if the 
eriterion 13 were Inhecessal from the vu ypoint 
of estimating an upper bound of error in the solution 
il would still bye a desirable conadttiol lo TLL Pose hh 
order to insure that the computed values differences 
with reasonable eas We shall further require that 
the term ol / InvOoLVIng the lowest ower of / 
shall approximat the magnitude of 7 to within a 
factor of two. The fact that a restriction is thereby 
imposed on A must be clearly kept in mind. In 
general a value of A small enough to satisf) 13) will 
not necessarily make / small enough to meet 


requirements for a given upper bound of error in 


I > se i le of & 





the olution In the instance when it does not, the the truncation term In other words. we shall 


problem posed to choose \ and / idiciously, so pare various choices of \ and f/, for a constant 7 
as to bring the truncating error within permissibl can be readily shown that for equations belong 
bounds Ww th the least amount of worl But if the type I, no other choice of h and AX keeping Z cons ; 
required accuracy ts rather low, it may well happen will be as good In this sense it is correct to ' 
that after / chosen small enough to satisfy 1 that X 1 is the best choice for formula (4 a 
the error measured by / ma niread be small ever, om Important point has been overlooked 
enough to meet all requirement In that case the have seen that A is not completely free, for kh m 
largest \ within the admissible range will, of course small enough to meet the minimum condition ' 
lead to the least amount of worl posed by (13 It may happe n that after A has | 
With these observation ve shall now attempt to taken small enough to insure that (13) is sat : 
tucdly thre cdependance of the truncation term on A the conditions for an upper bound ol error ! 
and on / Ideally, it would be di rable to classify solution can be satisfied with some range ol 
differential equations into several type recording to In that case, we would certainly ret a more a 
the value of } which appropriate tor equations solution by choosing X , but that would be 
belongime to the typ mn questlo \ compl tt accuracy that required and work could be save 
classification difficult to set down, but an attempt choosing a larger X\ We conclude that ever 
in this direction made DY Cor aeril two type equations ol Lv pe |. the choice of X ; will be 
Tyre Ll. Th [Vv pe characterized by the follow only if a relatively high accuracy is required ! 
ing condition solution Moreover, in practice the difference 


14 u(r,t) can usually be judged only from the in 
() 14 riven profile, at the time when A and /é are cho 
[ so that a safe interval A (rather than the maxu 
permissible for the given profile) is often choser 


T ole | ‘) le | J ] 
When +) hold it Toltow aa it may happen that the value of A considered ni 
sary may, as stated, bring the error withn 
au eu , ; tolerance limits for all choices of X 
Of Ju )/ ) Still considering svstems that be long Lo type | 
us now examine approximations of order four | 
( rT at to vst s belonging rt 
It ean hown that for o longing to type | expression for the truncation term 7, iS givel 
the terms of 7 } ft, do not involve f(z7./,%) orits 
m7" 9a) and the term of order four is 
derivative Thre bast homo enou equation be longs 
to this Ly yp 
Tyee Il This type | characterized by the con h*Ay(a,t,d }, ( A A j = AOvMlOnr 0 
dition that successive derivatives of f) with re 6 1? ° 90 Ou O 
spect to? are known to be much smaller than corre 
sponding derivatives with respect to usually of It can be readily verified that \7(\) is positive fo 
twice the order) with which they are associated in | ¢hoices of \. hence the leading term of 7, cam 
' a - as (7b) and (9b 5 
truncation term uch a a and , It quite be eliminated completely “us in the case of 
easy to find svstems belonging to this tvpe the 
: simpler approximation of order two. However 
numerical illustration given in section 5 belongs to 
type II may seek that value of \ which will make the w 
e 
Consider differential systems of type I, and let us work, Z=1/)\h*, a minimum, subject to a g 
start with the case when formula (4) is used. The | Permissible upper bound of | 7; By hypoth 
corresponding truncation term, 7 . 8 given by h is small enough so that A*A,(z,t,r approxim 
7a It has been observed by Milne, as well as by | the magnitude of | 7;.,,.,|; henee if Cis the permiss 
Thomas and others, that if \ is chosen as }, then the | tolerance of | 7, . we wish to satisfy 
leading term of 7; ,,, Vanishes, so that the order of 
Ty » » becomes four. With \=} and A small enough ou 
to satisfy the conditions for an upper bound of h*)| Ag(a,t,d hi M(x Ds <( 
. . . J 
] enough to meet the requirements for a given 


upper bound of error in the solution—a certain 
since (|O°u/Oz' Is independent of the choice 


parameters A and \, the inequality expressed in 
can be satisfied by taking 


amount of work will be done, which we shall measure 
in units of Z—1/)A*, as it has been shown that for 
a rectangular lattice the work is proportional to Z, 


approximately. If another \ and Ah were chosen, 
then we shall ask whether, for the same amount of htMnr C: <0 0 ~ 
work, and using the same formula (4), greater ac- Ou 


curacy can be obtained in the solution, assuming 
that accuracy to be measured by the magnitude of | Thus we seek to determine / and A, which minimi 
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to the condition 


WA) = ¢ 17 


17 / ( VA this value 
toZ and dilferentiat 1 
th respect toOA, We obt 


| / 


Substituting 
the resulting expression 
un the following condi 


ror a mn oO 


F—3M'(\)—4 MQ) 0, 


H 
' HOA i, S i) IS 
positive root ol ‘ \ O PHOS and 
rified that / \ positive, so that Z ts 
damimimum for this value of X Since X le 
n the admussible range t can be used for the 
itio, In conjunction with a suitable value of / 
ho Satisites LS nh practice 10.25 will 
al be used, becau i rational value of A Is 


enient 





vill bye mstruct ‘ Oo eNXamine thre followime 
ile. which gives AJ(A), A. and the work unit 7 
‘ fant value hi ALTO and various values of \ 
n the admissible range, in units of the corr 
ng quantities when J 
' y. an J 
1) (HOTS ] yticol  ,S2 
ool l«aisl $44 
OOO6o4 2 OOO} 2) 
Ow 1 visi 207 
OOLSS1 ] »HS/ +0 
003209 lL. 3o4/ rf 
O1l111 j 
O249U630 =. 7S3/ ] ro 
( ipared with a unit of work Z, for the case \= 4 
17, is required When A=}. The largest admis 
value of X 1s the poorest of all. from the view 


Again, a 
For the same magcnit uc 
and different 


of the 
rd of caution ts required 
the error in the leading term of 7; 


amount of work required 


hoices of X. the above schedule shows that with 
\ : hk can be chosen twice as large as that required 
hen A 4 If A is made twice as large, the eighth 


ference in the #-direction ts multiplied by about 2° 


d the differences must be reexamined to see whethet 
his larger value of the eighth difference still satisfies 

fundamental conditions imposed by (13) namely 
hat successive differences in the r-clirection bevond 
he fourth be numerically smaller than preceding 
nes Moreover, it has been pointed out. before 
that the maximum / corresponding to which (13 
satisfied may already be such that i i is within 
the required tolerance for all admissible values of \ 
In that case it is of course best to choose 3 or one 
to it 


Close 








The arises nee the 
approximation (4) has a truncating 
four when A and 7 . 

there any gain in using the approximation formula of 
higher ordet! The answer to this question ts com 
plicated by the fact that the 
ate a profile corresponding to the higher approxima 
mav be considerably longer than that for the 
simpler approximation. Much will depend not only 
on the computing instrument which will be used, but 
conditions 


following simplet 


que stion 


error of ordet 


also ol order tour is 


time required to venet! 


tion 


also on the complexity of the boundary 


If the IBM ecard programmed calculator is to be 
used and the boundary conditions ure not too com 
plex, the simpler approximation (4) can perbaps be 


renerated in only two-thirds of the time per profil 
more claborat 
Ther ure 
the highes 


tion worth considering. Let us assume we are dealin 


compared with thy fourth order ap 


proximation given in (5 however. com 


pensating factors which make uUpproXita 


with a case where the “best values of X are used in 
the approximation of order two and the one of order 


kor the 


of lattices points is inverse ly proportional toOA 


number 
ae hee 


two-thirds the number of lattie 
This 


hay 


oul same A used in both cases, the 


we shall use only 
points when the higher approximation is used 
compensate for the longer tims 
There is 
accuracy When the higher approximation 
coefficient of h*(0°u/du 
responding to the simpler approximation, is 1/540 
when A On the other hand, the coefficient of the 
corresponding term in the higher approximation with 
nN ; is only 0.000694; or less than three-eighths that 
of the simpler approximation Furthermore, in 
cases Where the maximum permissible /& is such that 


would abo il 


take to generate each profil however, a 
vain in 


used For the In (Ua col 


T, is already within the required tolerance limit 
for all values of A, we may be able to take A= 2, o1 
when the higher approximation 
Such a choice of \ would cut the number of lattice 
+ that required for the simpler approxima 
tion, with A=4, and that might mors 

sate for the vreatel ciffic ulty In venerating rt) by 
the higher approximation 


close to if Is used 


points lo 


than compen 


that be- 


deriva 


Let us now consider differential svstems 


long to type I] ses by hypothesis thr 
tives in the ¢-direction are negligible compared with 
leading 


those in the x-direction 7b) shows that the 


term in 7 is practically independent of A, and 
choice ol \ 
it is reasonable to take the maximum admis 


Similarly, 


this term cannot be eliminated by any 
Henee 
sible A}, or on 
term of 7, 


no optimum A 


close to it the leading 


is complex in structure, and there is 
that 


functions falling under this type 


stands out as suitable for all 
For a given prob 
lem, it may be possible to compute estimates of thi 
various terms that contribute to the truncation error 
Or if the problem involves a family of parameters 
in the boundary condition, then available results for 
some members of the family may lead to an optimum 
Nand A for the 


choice of remaining ones 














high I 


For 


upproXimation 


falling 


formula usu: 


equations nae! 


Lv pe 


illvy has distinct advan 


tages over the simpler one of order two For in this 
case the truncation term is necessaril of a lower 
order of magnitude in the high: ipproxXimation 
formula than in the simpler om moreove a large! 
mesh ratio, namel can be ta ! Whenever the 
boundary conditions are such that the coding prob 


lem Is manage the 


commended 


able gher approximation formula 


to De re 


It might be well to remark that in some cases the 
olution mav be an oscillating function of Dut not 
ol f Such functions may fall under tvpe Il, and a 
Fourter approximation may actua be better than 
a finite difference approximation However there 
sno reason to expect that all, or ever rood portion 
of functions belonging to tvpe Il, can be most simply 
treated by Fourie approximations \ finite differ 
ence approximation such as (4) or (5) is often pre 
ferred, because of its simplicity, to othe types ol 
approximations (by Fourier series or perhaps “im 
piieit difference approximations us concern 
here, as stated before, is to study the dependence of 
the solution on the choice of \, afte thet t) or 9 


has been selected as the formula 


approximatlor 


4. Method of Improving the Solution from 
Two Difference Approximations 


he method to be expla ned below has been 
presented in [1 by Hartree and Womersley, who 
ascribe the idea to L. F. Richardson [5 In 1] the 


method s applied to a somewhat different type of 


approximation a” mixed difference differential 
scheme, suitable for computation by a differential 
analyze The results will here in extended to 
difference approximations of anv othe 

Let us suppose that values of have been gen- 
erated by an approximation torm ila corresponding 
to a true solution and let us assume that it ts 
possible to writ 

h’¢ fxr / ( f x 19 

where the functions ¢ fr) are independent of & 
It can be shown that corre sponding to a wide range 
ol boundary, conditions the expression [or does 
assume the form 19) for both approximation 
formulas (4) and (5 

Consider now the case when has been gen 
erated by the use of an interval A, in the «-direction 
and a mesh-ratio \ Let these values of be 
designated by h Now let another computa- 
tion be made, based on the same mesh ratio \. but 


at an interval A, in z, where h,= ph,, 0< p< 1 These 
values will be designated h It is presumed 


that the same approximation formula will be used in 


by 


‘} ure brack licate the literature refer ! this pape 
I bsequ sssions we shall write eh r o(hi.hs) t ndicate 
h oT Pm.al(h,A when no ambiguity like t ‘ } 
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both cases Clearly 

h hi l / i 
By hv pothesis we have from (19 
/ } h’¢ / ( 


/ ( 
Similarly, remembering that / ph, and that ¢ 
independent ot / 
f p hit | ( 
ia | 
Niult pl 70) by p” and subtract from (2] 


ransposing some terms and divid 


W here 


ss defined 


the proce 


We shall refer to in (2 is 
“p * correction In 1} itis recommended that 
taken as 3. This ts a desirable choice for num 
work, since every other value in the z-direction at 
smaller interval is available at the larger inte: 
Similarly, every fourth value in the ¢-direction 


be available at both intervals There many wi 
of applving corrections of this type One 
exampl to generate four values in the 
at the finer interval, then generate the 
values at the larger interval; apply the p 
to the last profile and use these corrected \ 
initial data for generating the next profile. Sucl 
use of the correction scheme would, of course 
h h for ot! 

for those po 
The coding 


are 
wa 
Is f-clirect 
correspond 
correct 
alues 
requ 
interpolating for values of every 
can correct only 
both intervals 


value of 2 we 
that are available at 
this method would be complicated. The simples 
wav of using the correction to actua 

make two completely separate computations for 

required values of n, and then to apply the corr 
tion process only to those functional values that a 
Often the last profile generated is of mos 


since 


st heme 


required 
interest, and in that case it may be enough to corr 
the values on the last computed profile only 

The process furnishes a powerful check on the co 
vergence rate of the approximations, and sinc¢ 
work done at the larger interval is one-eighth th 
done at the smaller interval if p= 4, the added lab 
when the two computations 
Moreover, the coding is the sar 


not too costly , 


is 


carried separately 

















schemes For the specia 


oth approximatiot 


when 2, (22) and (2 reduce to 


WiC 24 


nn ‘ : | 
happiied in the numerical 





yes ven in se ( vith highly satisfactory 
s. th Hartree d Womersley give sufficient 
ons ol tha natu ol the boundary lol the 
“ o be Vali l method can probably bye 
OV i widel ss of tunctions than those 
ad in 
s sufficient } will alwavs be a 
nent ovel the truncation term 
order oO ide ou h, J Th 
) irises LO ! st A be to insure that 
shia be ! rove vent ovel if Let 
hneacol | / how by kien! ed by 
vhen assor th an interval A in the 
Ons vrhe 
hj j 


e sSatisiied \ nume cally smallet runcation 
vill usually | - d with a smaller total 
Kor the approx ition formulas considered 

: truncation term ts of the form 


larantees 


no similar statement was made about 7 h 


{ 


irom 


































where 


shown by summing 


/ h, | that 


then from (27) and 0) it can be 


the absolute values of the 


27) is satisned When ; and 7 2, the condition 
29) implies that p If ; + and +, then p z> 
or What Is equivalent, successive terms ol / at 
the larger interval are required to be just smallet 
numerically than preceding terms. Under such con 


ditions v(h, he) will always be an improvement over 
j j 

In the foregoing t has been assumed that can 
be computed exactly by the prescribed formula, and 
that all initial and bounda conditions are exact 
This can seldom be realized in practice and there 
vill be rounding errors committed at every step of 
the computation, due to carrving a fixed number of 
decimals or significant figures in the computations 


he cumulative effect of such errors can perhaps bi 


studied statisticall or upper bounds tor errors of 
this t pe |Z) « in be found In then merical examples 
viven in section 5, the cumulative round-off error 
vas vi ~ ll, af 00 steps u 


5. Numerical Examples 


em selected for analvsis was the following 


t? j 
0 
This pa ‘ ir form was chosen because t repre 
sents a cause Nn which aerivatives W th respect to f 


smaller numerically 


aerivatives with 


are much than corresponding 


; 


respect LO JZ with they are 
iated in truncation terms such as those of (7b 


Il of section 


It is known that the exact solution to 1 


which 
nssor 


The t 


like re ntial SVstem bye longs to tvp 


choice ol B® noniinedal form a the dulteren ul 


The 
‘ 


equation was deliberat It was made in order to 


study the where a 


cumulative 


error Ih Cases 


siderable number of operations, that involve the 


approximate values of u, have to be performed at 


en h Ste p 







The following five sets of solutions were generated 








Parameters / 0.05 / 0.00125 ? 

Range of f: O0<?t<0.12/ 

Number of lattice points 20) OO ? OOO 

Initial values and the computed values of v,., for the last profile are given in table 1] 

| \ } 
0 ) .f | L1.¢ 
t} 
O.05 00125 
y LO. 0.04 
\ 0.125 
0. 00 0 SO00000 0. S720400 0. 8746427 0026027 
OS 8155520 SS6H4322 SASSRO7 OO2Z4575 
10 SOHOSTLS6 Q285124 9305912 OO02Z07S88 
15 9318664 9950813 Q966S819 OO L6006 
20) 1. 0237187 1. OS1L9OL0O5 L. OS30512 OOLL407 
25 1316470 1. 1848335 lL. [S55005 OOOTSTO 
4) 1 2aLT OR; lL. 3003062 1 3007690 OOO4628 
35 1. 3812585 1. 4255420 1. 4257892 0002472 
1 1. 5178933 1. 5584295 1. 5585249 0000954 
5 1. 6601581 1. 6973964 1. 6973876 QOOOOSS 
a0 1. SO69311 1. 8412718 1. 8411952 00007 66 
55 1. 9573898 1. OROLSILO lL. YS90638 OOOLLS1 
60 2 1109240 2 1404665 2. 1403270 OOOL3S95 
65 2. 2670741 2. 2946235 2. 29044770 OOOL465 
70 2. 4254896 2. 4512642 2. 4511221 0001421 
75 2 5858903 2. 6100866 2. 6099568 0001298 
gO 2 7480902 2. 7708507 2. 7707399 OOOLLOS 
85 2. 9118937 2. 9333664 2. 9332788 0000876 
o0 3. 0771740 3. 0974789 3. 0974182 OOO0607 
95 3. 2438210 3. 2630627 3. 2630315 0000312 
1. 00 3. 4117444 3. 4300145 3. 4300145 0 


EXAMPLE 2 


Formula: The same as in example | 
Parameters: A=0.1=h,; k=0.005; A=} 
Range of t: The same as in example | 
Number of lattice points: 10 25=—250 

Initial values and the computed values of v,,,, for the last profile are given in table 2. 
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} 0) 

i 

‘ ) 

> 

s 2 

, > 
PI 


Formula 
Paramet« 


Range o 


Number of lattice 
Initial values an 


{ 


formula for examples 


| 


i 





0. SOSSS7S 
9221349 
oO, 84844 
2080454 

S1L507 
S415191 
LOSUS4 

2 4517066 

2. 7411921 
OOTH676 


§ 4300145 


i example | 


/ 1/1176; 


ane as in example | 
pomts 14 147 


computed values of 


to 4 comes trom 


SOMO) 
8314955 
Y203245 
OS531018 
2164087 


LOOSSO0 
5OSSTS 1 
8069311 
0228433 


2. 2446210 


2,4711277 


» FOLSTSS 


2. 9354176 
3. 1722397 


3. 4117444 


O.005 


305912 


(4042; 


OS30512 





PO5S 


1 


;O076900 


9585249 


for the 


SHO2554 
SOS7T5H40 
YS246030 
LOST 660 


2H56385 


1438567 
6371338 
S413398 
0537601 

2725815 
1965649 
248430 
O567941 

1919624 
100145 








OOS4156035 
OO4 5668 
OOTS2Z56 


0005, 42 


OOOS234 
OOOS669 
OOOSS45 
OO04522 


OOO24904 


8746428 
9035402 
QR5S00 1 
LLOS660 


2667479 


1443213 
6372086 
$411953 
S35080 





(22948 


1962917 
7246151 
9566326 
I9ISSI1 
1300145 


() Sy, 4, 


Y2004S 
OS3052 
S00, OOD 
»OSH2Z2 


| SOR 


14032 


oOo741¢ 


viven in table 


OOSSS74 
OO4TS53 
0034361 
0021000 
OOTLOO4 


OOOU4646 
OOOOTAS 
OOOL445 
OO002521 


OOO02867 


OOO2T 32 
0002279 
OOOT6LS 
OOOOS 13 


boll lGs 


> 7.01369 





OOOLL4, 
(M04, O 
OOOOOLS 
OOS) 


OOOO 5 


OQOO0059 
OQOOOO30 
OOOO054 
OOOOOS0 
OOOO022 








EXAMPLE 4 
Formula 


It Will be convenient 


form 


\ a 
where 
a a MN: oa 


Parameters: A=—0.0625 


} 


Range of ¢ O<f<0.12S8SS 
Number of lattices pots 
Formula 2) 18 a modilicat 


tamed by dropping thi terms oO 


However. the second deriva 


l 





OOOO 
OO25 
1250 


1S 


3750 
1375 
SOO) 


5625 


7500 
S125 
8750 
9375 
OO00 


f 


nvolving p>2 


| 172 


, 
It 


»A 
3584 
is ob 

) 


with respect to 7 1s 


SOM MMMM) 
8242006 
SU34221 
QO900767 


r3ih4d,0 


2PR33S862 
14S, Si 2 
6241314 
SOHUR 11 


9955052 


LSS8OH06 1 
S856893 
58589093 
7S88O57 
Q9043567 


2020364 


1117444 


table l 


0.0625: 


8767367 
SORSOZH 
9627974 
0615954 
1872134 


3326403 
4925905 
6633207 
8422337 


027: 





2179277 
$125343 
6106854 
SILIS924 
OLST7S64 


2220847 
1305684 


approximated by differences, including the fi 
order. The resulting formula is almost as ac 


as (5) for this problem, since derivatives with r 
to ¢t are very small numerically In view of th 
that (5) was modified, it is necessary to examin 
to determine the admissible range of A It is 
from (33) that all the coefficients a, cannot bs 
positive by any choice of A; hence the stabil 
teria given earlier do not apply However 
been shown in [2] that if there exists a positive 
ber ./, independent of « and f, such that the ec 


cients da, of (33) satustiy 


» > a, exp qy)|\S exp (= V/ J"), lor y 


then the bast homogeneous differential equatio 
stable From that result we may then deduc 


stability of 2) or (33 


, qua , 
yh 
, 
4) ‘ j | 
0.0014: \ 0.3584 
At 0.1288 
0. 8768124 OQOO0OT 57 
SOS0475 QOO0545 
96828100 OOOOT25 
1. OBLSS04 OOOOLSO 
1 1871920 OOOO2 14 
1. 3326211 OOOO 192 
1. 4925762 OoOoo1L43 
1. 6633105 OOOO 102 
1. 8422269 OQOO006S8 
2. 0275206 OOOOOAS 
2. 2179247 OO00030 
2. 4125324 OOOOO19 
2. 6106848 OOOO006 
2. 8118925 OO0000 I 
3. OL57871 OQOOO007 
3. 2220858 OQOOOOT 1 
3. 4305684 0 
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MPLI 
Formula: ‘The same as in example 4 
Parameters h 0.125: k—0.0056: \—0.3584 
Range off: Thesameasinexample 4. Results ari 
PaBLe 5 N, j 
’arame } 0.125 
/ ; ; } 
h give ible 4 
t) 
1 SOOO) Oo. S75 sH2S i) S7OS124 
Ld 25 8934221 9623470 9628099 
) 1. 1316470 l. 1872934 1. 1871920 
1. 4487872 1926595 1. 4925762 
i} 1. SOOUS11 §$4222406 lL, 8422269 
2 2. 1 SS86961 2. 2178929 2. 2179247 
i) 2. SS5S993 2. 6106245 2? 6106847 
R75 2 9943567 § OL5S7128 3. OLS7TS869 
O00 §. 4117444 3. 4305684 §. 4505684 
In all the above examples u and u 
cenerated Fourth differences were cenerated 
ll the examples, even where the computing for 
a did not call for them. These fourth differences 
jotted down bv the operator of the card pro- 
ammed IBM calculator, and anv lack of continuity 
the diiferences was a warning that the machine 
| 





not Turn tioning prope ris The ope rations were 
a board wired to perform eight- 


mi out with 


calet 


multiplication Phi lations were carried 
the fullest possibl accuracy, that is, to seven 
ials in and to eight decimals in some of the 
sidiary computations for This 


iracyv is far in excess of the truncation error, for 
ull values of 


\ ‘‘p 
available profile 


’ correction W is applied to the values of the 
as explained in section 4, based 
] and 2 The 

Similarly, a p* correction 
profile in examples 4 and 5; thi 


Sarre given in ¢ xample ) 


the entries in examples results are 
en in example 2 


1) ed Lo the last 


; 


Was 


5.1. Observations 


It is to be noted that in spite of the fact that in the 
profile differs from the true values by 0.0026 
example 1 and by more than 0.01 in example 2 
values of r h he resulting from the p correction 


to within 0.00011 or better It can be 


correct 








( , hed 
7 
0.0056 } 0.3084 
‘ j , 
p I f 
\ 0.1288 
hyd h,} 
0014496 0). S7TOS282 OQOOOTSIS 


9628274 
1872081 


OOOOLTS 
OOOOT61 
(MO00907 


0004629 
OOOTOL4 
OOOOSS3 


1925859 


QOO00027 1. 8422340 OOO0O07 1 

OOOOS TS 2 2179300 OOOO053 

QOOOHO2 2 6106895 OOOOO4S 

OQOOOT 41 , OLST913 OOOOO4 4 
t) ‘) 


verified from the tabulated entries that 


~he th 


It follows that to fain an accuracy comparabl lo 


that of vr(h,.A.) without the ‘‘p*” correction, an in 
terval of / 0.01 would be needed, and hence the 
amount of work would be 125 times that used in 
example | 

The upprovement due to the p eorrection in 
examples tf and 5 is not quite so striking However. 
even here there is considerable improvement for 
small values of and it must be remembered that 


an A* term in the trun 


eliminated by the p' cor 


the formula used does include 


cation error, which is not 
rection, although the error from this term is somewhat 
than v(h in 


than v(hy.h 


However. for large! 


value 


lessened 


example 4 


Is closet to the true 


An explanation for this may come from the following 
onsiderations (ssume 
then it follows that 
" p h ( p | p 
u hyh he 
l p | f 
o(1 ' 
r! h5¢ 6 
] p 





A? and h 
than in 7) How- 


noticed that the terms Involving 


smaller in > 


It is to be 


are rmhumerict 


ever, the term involving /° is somewhat larger, if 
} and all subsequent terms will be larger In 
fact, it can be readily verified that if a ‘ * correc- 


tion is applied, then compared with a term involving 
is the corresponding term 


/ ( in ? there 
{ l { | p il ( ih | hy Jf which is, 

isually the leading 
aller than the term 


approximation 


ot CoOcdrsé 


after the 


numerically larger Du 
term elimination is sm 
eliminated, namely, A’ p’¢ The 


ais applied inh examples tf and 5, is rather unusual in 


that a verv small term of order two is left If all the 
terms ol o>) eae small numerically it may happen 
that the combination of the leading terms in u 

hy he mav be somewhat larger than in wu h 
In such cases, however, the difference between r(h,.A 
and h.) will itself be small: hence, although it may 
not be known which is the better answer, it is to be 
expected that the order of magnitude of the error in 

hd is not greater than / hy fh 

The last place ot u f)is not cuaranteed he nee we 
can judge the rounding error only to the extent that 
the sixth de« imal place is affected It is to be noticed 
that in « xam pli t, which seems to be the most ac- 
curate, the difference u,, , is systematic as far 
as sign is concerned This is evidence of the fact 
that no lara rounding off error accum ilated, after 


differential 
numbe r of 


even though a nonlinear 


and a 


Loo steps in 7 


equation was used considerabl 


arithme tie ope rations were pe rformed at each step 


] 0. YLOB6SO2 

0 SOO0000 2370047 
] 9203245 1757597 
2 1. 2164087 OSG60852 
3 1. 59SSTS1 0420958 
! 2. 0228433 0240192 
a) 2. 4711277 


The fourth differences are not numerically smaller 
than the third differences, the first entry in the 
fifth-difference column is very much larger numeri- 
cally than the second one in the same column. Such 


and 


a pattern is a warning that the interval is too large. 


r r t? - w 
ia he Eftect ; na a Hicher F Arnnroxyimation F 


Let us compare the results in examples 1 and 4 
The coding for example 4 is somewhat more compli- 
cated than that for example 1, but since only 1,472 
lattice points were used in example 4 compared with 





eXalp 


Let us compare the error pattern in 
where 


where A=4, with that of example 3 
In spite of the fact that somewhat more wor! 
performed in generating values in exampl 


results are not as eood although the error in 
examples is in the same decimal place; but the 
in example } is about twice as large This, in fa 


precisely what was to be exper ted For as obss 


in section 3, in cases coming under type I] 
derivatives with respect tof are relatively smal 
error 18 not appreciably aifected by XA: henes 
The ratio of the two va 


almost proportional toh 
numerical illust: 


of A? is (20/14)?~2.04 The 

verifies the observation that A ‘ is not necess: 
the best mesh ratio for the formula given in 
The boundary conditions and the form of 


determine the type that the diiferential svstem 
longs to, and it is only after the problem is stu 
from this viewpoint that a suitable choice of ) 
be made 

In connection with example 3 
X 2 it was not desirable to generate values at do 
the interval in the z-direction, for purposes of ay 
ing a ‘‘p*”’ For if the formula giver 

} were used to compute we would vet a 

tive Even if the 
were 2 knowledge of the more acc 


eorre spond ne 


eorrection 
value of true ans 
known 
, from the smaller value of A would 
incorre 


not 

values of 

negative value of is 

instructive to write down the differ: 

doubl 

and to use a valu 
The values of 


us that a 
It will be 
of 3 the known 

interval used in example 3 
computed from formula (4 
some of the ciiferences are given below 


initial profile) at 


0612450 
OS896745 
0439894 
OLSO766 


0284295 
O456851 
0259128 


2,000 in example 1, the over-all amount of wo 
about the same in both cases. The result after 
steps shows that the higher approximation form 
gives very much better results. To secure a ma 
mum error of 0.00008 in r,, , with the approximat 
used in example 1, it would have been necessary) 
use an interval A of about 0.009; hence 170 times 
amount of work would have been necessary. Hi: 
ever, if the p’ correction were applied, the compa 
tive results would not be quite so unfavorable to 
Assuming that a p 


col 


simpler approximation 
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were appli d to two computations by the simple1 

roximation, and a p* correction to results of th 
fied fourth-order approximation, the latte! at a 

ld still five significantly better results om ter me 

tional decimal place in fact To secure com- 

ible accuracy by the simpler approximation, it | 7.2. Relations Between Derivatives 

ld be necessary to multiply the amount of work r-Directions 

the factor 10 The conclusion is inescapable 

B i i 1 
the higher approximation 1s worth while in 

s where the boundary conditions do not introduces 

ilarities in the higher derivatives, and when the 

ng probl m is Manageabl 

Phe author acknowledges eratefully the many 

structive suggestions which were given by D1 if / 

tz John during the progress of the study 
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The third term on the right-hand side of (49) wa 
dropped in example 4, since its magnitude would 
affected only the fifth decimal place at any point 
term could have been obtained, if desired, by using 
and (37), and then setting m l and » 2 in (49), t 


three linear equations for the unknowns ‘ and h 


J 


7.4. The Boundary Condition at the Termi 
Points (z,, ¢) 


The difference equation defined in (4) needs no 
treatment for the boundary conditions at z,, where 


since (s l.n) is the last lattice point at which p,, , 
erated, and &2 » is fully defined However, whe 
fourth order approximation is used, as in (5), an expre 
is required for v,.; ,, or for ®u, Since u(zrz,lt Bit 
all its derivatives in the ¢t-direction are a imed to be k 
the differential eq (1) can be used to obtain the 
partial derivative in the z-direction, in terms of Ou 
f(xr,u The relations between derivatives and difference 
then be used to obtain &u,, In examples 4 and 5 
fourth difference in the z-direction was essentially 


Le 
Z 


r,; hence the last fourth difference was replaced by & 


Los ANGELES, CALIFORNIA, September 18, 195 





Thermal Properties of Some Butadiene-Styrene 
Copolymers’ 


George T. Furukawa, Robert E. McCoskey, and Gerard J. King 


he thermal properties of 41 


and 122° F butadiene-st yrene 


copolymers containing 8.58 


percent of bound-styrene were investigated by means of an adiabatic calorimeter from 16 


to 330° K The 41° F 
210° to 285° K The 122° F 
transformation temperatures of the 41 
respectively 
have been 


pute heat 


matior temperature 


were med te con capac 


intervals 


l. Introduction 


The literature is copious with studies involving 
physical and chemical properties of various rub- 
r polymers prepared from different recipes. These 
vestigations have principally in mind the correla- 
m of the observed properties with the structure of 
the polymers—the structure or the composition being 
termined by chemical methods, derived judiciously 
om the units that compose the polymers, or deduced 
from the various observed properties themselves 
The situation is complicated by many variables in 
the recipes, such as polymerization temperature, 
ercentage conversion, catalyst, modifier, and others, 
both controllable and noncontrollable, which affect 
the internal structure of the final product. For 
xample, in the polymerization of dienes many of the 
ibove factors influence the degree of 1-4, 1-2, cis, 
trans addition as well as cross-linking. Further- 
lore, IN many instances the noncontrolled variables 
outweigh the controlled variables in determining the 
physical property 
Some of the phy sical properties of rubber poly mers 
have been known to be dependent to a certain extent 
ipon their previous thermal and mechanical history 
For example, in the heat-capacity measurements with 
Hyear O. R.-15 [1]? and GR-S [2], the direction of 
he temperature drift in the temperature range of 
lass transformation is shown to depend upon the 
rate at which the polymers were cooled prior to the 
easurements. The temperature drift observed can 
explained in terms of the long relaxation times 
ncountered for certain degrees of freedom in the 
polymeric substance to reach equilibrium. At the 
ower temperatures a given degree of freedom can 
be essentially frozen-in, and at the higher tempera- 
ires the attainment of equilibrium for that degre: 
of freedom can be very rapid. The temperature drift 
s observed in the intermediate temperature range 
when the relaxation time for the degree of freedom is 
of the same order of magnitude as the time for the 
the temperature 


heat-capacity measurements. In 


copolymer was found to crystallize ir 
copolymer did not exhibit any 
122 } 
rhe effects of heat treatment upon the heat capacity and the glass-transfor 
studied he 


ity, enthalpy 


and 


the temperature range fron 
ass- 


crystallization The g 
copolymers were 200° and 193° k, 
measurements 


5-deg 


the heat-capacity 
to 330° Kk at 


results of 
and entropy from 0 


range of glass transformation the relaxation time for 
certain degrees of freedom in polymers comes within 
this magnitude. Similarly, as is pointed out later, 
the crystallization and melting processes in polymers 
have a relaxation time within this range. Thus, when 
the polymer is cooled rapidly through the glass- 
transformation interval (or through the erystalliza- 
tion interval) a certain portion of the polymer mole- 
cules fails to undergo transition from the higher- to 
the lower-energy states. During the heat-capacity 
experiments and when heating to tempecatures near 
the glass-transformation interval, the frozen-in de- 
grees of freedom become gradually excited and transi- 
tions occur within the polymer in the direction of 
equilibrium. In this case, the polymer molecules 
tcansform into the lower-energy states with the conse- 
quent liberation of energy and upward temperature 
drift. On the other hand, when the polymer is 
cooled slowly through the glass-transformation in- 
terval, more polymer molecules have a chance to 
transform into the lower-energy states. Upon heat- 
ing to the glass-transformation range, the polymer 
molecules transform into the higher-energy states ac- 
companied by an absorption of energy and a down- 
ward temperature drift 

The glass transformation, which is associated with 
many polymers, supercooled liquids, and other non- 
crystalline solids, arises from the freezing-in or the in- 
creasing length of relaxation time for certain degrees of 
freedom and is characterized in calorimetry by a steep 
decrease in the heat capacity with decrease in tem- 
perature. The theemal coefficient of expansion de- 
creases rapidly in a similar manner, and there are 
many other physical properties that are affected simi- 
larly within the same temperature range. The glass 
transformation is not as sharp as first-order transt- 
tions. but 10-deg interval 
As the relaxation time for equilibrium in this temper- 
rder of magnitude as the 


Is spread over a 5- to 


ature range is of the same 
experimental measurements, the heat-capacity values 
obtained are dependent upon the rate of cooling 
through this interval as well as the time awaited for 
temperature measurements. Consequently, the heat- 
values are somewhat scattered 


eapacity 
other natural and 


Natural 


rubber and svnthetie 
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polymers are known to crystallize over a wide range 
of temperature There is only a short-range order 
within any crystalline polymeric substance and along 
a single polymer molecule both crystalline and amor- 
phous structures exist simultaneousl\ Recent X- 
ray studies [3] have shown that the crystallinity in 
natural rubber increased with storage, and that even 
after 30 vears of storage the crvstallization process 
had not reached the state of equilibrium The ervys- 
tallization and melting processes of polymers are ob- 
served as temperature drifts in calorimetry. The 
lower temperature limit of this crystallization range 
is set by the lack of thermal motion to attain crystal 
configuration; and beyond the upper limit, the ther- 
mal motion is sufficient to make the time spent in a 
crystalline state short and the polymer is amorphous 
Thus, the degree of crystallinity of the polymer at 
temperatures below the crystallization range, is deter- 
mined to a certain extent by the rate of cooling 
through this range 

This report deals with the heat-capacity investi- 
gation of two butadiene-styrene copolymers contain- 
ing 8.58 percent of bound-styrene and prepared at 
41° F (5° C) and 122° F (50° C). These two mate- 
rials were investigated to determine the effect of the 
polymerization temperature upon their relative crys- 
tallizability, glass-transformation temperature, and 
heat capacity. The effects of heat treatment upon 
the heat capacity and the glass-transformation tem- 
perature have been studied, and the results of the 
measurements were used to compute heat capacity, 
enthalpy, and entropy of the polymers at 5-deg 
intervals from 0° to 330° K 

The calorimetric method, as applied in this inves- 
tigation, is sensitive to changes in the polymer in 
which thermal energy as little as 0.03 abs j min™! 
is evolved or This method was used to 
detect temperatures at which thermal effects occurred 
and at which transitions took place. The tempera- 
ture range of crystallization in the 41° F copolymer 
and the glass-transformation temperature of both 
copoly mers were determined 


absorbed 


2. Apparatus and Method 


The details of the calorimetric apparatus and 
method used in this investigation can be found in 
the previous reports dealing with 1,3-butadiene [4] 
and diphenyl ether [5]. A small quantity of helium 
gas, sealed in with each polymer sample, supple- 
mented the vanes in enhancing the attainment of 
thermal equilibrium. The temperature measure- 
ments were based on the platinum resistance ther- 
mometer calibrated at this Bureau in accordance 
with the 1948 International Temperature Scale [6] 
and between 10° and 90° K on a provisional scale [7], 
consisting of a set of platinum resistance thermom- 
eters calibrated against a helium gas thermometer. 
The temperatures expressed in degrees Kelvin were 
obtained by adding 273.16° to the observed temper- 
atures in degrees Celsius 

In order to study in what way the heat treatment 
affects the heat-capacity results, two series of experi- 


ments 
from 80 


were carried out in the temperature 

to 280° K In one series, the polym«e t 
shock-cooled by immersing the calorimeter in | 
nitrogen with helium gas in the space surrouw 
the sample container By means of this proces 
the sample was cooled from room temperatur 
about 90° K in 30 min and to about 80° K 
additional 30 min In the second series of mea 
ments, the polymer was cooled slowly by main 
ing a high vacuum around the sample conta 
The cooling process was prolonged by SUCCeSS 
using dry ice and then liquid nitrogen. In s 
cases the 41° F polymer was left between 200 
230° K for a number of days to increase the ery 
linity by allowing a longer time for crystallizat 
After cooling the polymer to the desired temp 
ture, the heat-capacity measurements were Mad 
progressively higher temperatures, the final tempx 
ture of the first measurement being the initial t 
perature of the second measurement, and so fo 
up the temperature scale. As heat-capacity m: 
urements with simpler substances [5, 8] show 
normal temperature equilibrium was establishe« 
6 to 7 min after the heating period, any persist: 
temperature drift after the eighth minute was con 
sidered to arise from the thermal effects in the 
polymeric material. 

The net heat capacity, or the heat capacity of t! 
polymer, was obtained by subtracting the heat 
capacity of the empty container from the gross heat 
capacity. At the lower temperatures at which th 
slope of the heat-capacity curve changes rapidly, th: 
temperature Increase per heating period was made 
as low as 2 to 3 deg to make insignificant the curva- 
ture correction to the experimental heat capacity 
As the curvature does not change rapidly at the 
higher temperatures, larger temperature intervals of 
7 to 9 deg were used. The heat capacities tabulated 
at 5-deg intervals in tables 4 and 7 were obtained b 
graphical smoothing of the net heat capacities. T! 
heat capacities in the transition regions were obtained 
by methods later described 


3. Samples 


The two 90/10 butadiene-styrene copolymers 
vestigated were obtained from different 
The history of the rubber polymer, designated X—45 
and obtained through the Copolymer Corp., Bator 
Rouge, La., is not completely known. This copoly- 
mer was emulsion polymerized at 41° F (5° C) with 
the initial butadiene-styrene charge ratio at 90/10; 
the emulsifier was Dresinate 731 and the modifie! 
was probably Sulfole B-8. The reaction was “short 
stopped”’ at about 55-percent conversion, using 
ditertiary butylhydroquinone. No further informa- 
tion is known about this copolymer. The GL-658 
copolymer prepared at 122° F (50° C) was obtained 
through the University of Akron-Government Labo- 
ratories [9], and the details concerning this material 
are found in the reference given. The copolyn 
was emulsion polymerized with the same initial 
butadiene-styrene charge ratio as the 41° F copo 


sources 
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using, however, as emulsifier soap flakes (SF), 
as modifying agent normal-dodecyl mercaptan 
DM The polymerization reaction was activated 
potassium persulfate and ‘“‘short-stopped”’ at 
it 72-percent by means of hydro- 
one 
The polymers were purified by M. Tryon of the 
ber Section of the National Bureau of Standards, 


fa bl 
The 


conversion 


the purification procedure was as follows 
terial was first in benzene and centri- 
ed to eliminate solid particles. This solution was 
led dropwise into well-stirred methyl alcohol to 
gulate the polymer. The coagulated polymer 
; then redissolved in benzene. The solution and 
sequent precipitation processes were carried out 
e times. This procedure eliminated fatty acids 
their salts, stabilizer (phenyl betanaphthylamine 
low-molecular-weight polymers. After 
solving the polymer for the fourth time, about 
percent of phenylbetanaphthylamine was added 
the benzene solution as a stabilizer. Subsequently, 
benzene was removed by pumping on the frozen 
ymer solution held at the dry-ice temperature, 
iving the stabilizer behind with the poly mer. This 
Ss vielded a spongy mass which was pressed 
to ‘-inch sheets. The polymer samples were 
eived from the Rubber Section in this form, and 
vy were cut to size to fit between the vanes of the 
mple container. The samples were pumped at 
h vacuum for 3 days at 25° C and 1 day at 50° C 
eliminate volatile impurities such as moisture, air, 
benzene. Following this, the container 
led with a small quantity of helium g 
iss of the X-454 copolymer investigated 
¢, and that of the GL-658 was 40.349 ¢ 
\ portion of the purified polymers was analyzed 
The results are given in table 1 


dissolved 


some 


was 
The 


was 


gas 
0) 


oO 


composition 


The styrene contents given in the table were com- 


ited from the carbon-hydrogen ratios corrected for 
mercaptan content. The phenylbetanaphthyl- 
nine was added after the analysis 


PABLI 


4. Results 
4.1. X-454, 41° F Copolymer 


treat- 
and 


The experimental details regarding heat 


ents, temperature range ol the measurements, 


emperature drifts with the N-454 copolymer are 
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served Heot 


summarized in table 2. The observed heat-capacity 
values are given in table 3 and plotted in figure | 
The heat-capacity measurements In runs 1, 2, and 3 
were made after shock-cooling the polymer from room 
temperature to that of liquid nitrogen. In 
experiments the upward temperature drifts, observed 
in the range of temperature from about 170° to 280 
K, are attributed to slow changes in the polymer 
to of In run 2 
temperature drift at 220° K was much a 
deg per min 40 min after the end of the heating 
period As shown in figure 1, the heat-capacity 
results in run indicated by filled-in circles 
temperature interval from 

Then, there follows a sudden drop 
in the heat capacity The 
capacity rise the interval 195 205 
attributed glass transformation The 

transformation temperature is taken to be 200 
The sudden decrease that follows is caused 
crystallization of the polymer, for apparent 
capacity falls off with ervstallization The 
following this is from the melting of the polymer 
crystallites. The results (fig. 1) of all heat-capacity 
measurements are relatively high in the temperature 
interval 230 280° K, and during the heat 
capacity measurements upward temperature drifts 
observed in this It that 
irreversible melting took place in the vicinity of the 
central during the heating 
the temperature invariably higher 
the heater For this reason the heating 
the power, somewhat affects the heat 

During the equilibrating period some of the melted 
crvstallites reerystallize, and crystallites 
form with the increased thermal energy now avail- 
for orientation The apparent high heat ca- 
pacity indicative of a relatively larger amount 
of melting in the particular heating interval in 
comparison to crystallization. Thus the amount 
of the polymer melted is greater than the amount 
crystallized in the temperature range 230° to 280° K 
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although the observed upward temperature d 
indicate crystallization to be taking place. On 
other hand, in run 2, as pointed out earlier, th 
heat capacity just above the glass-transforma 
temperature arises from a large amount of eryst 
zation. As temperatures used in heat-capa 
calculations were determined during temperat 
drifts, the heat-capacity values are different { 
the values that would be obtained if thermal . 
librium had been reached The temperature e 
librium, however, was awaited sufficiently lon 
that the general nature of each heat- apacity \ 


TaBLe 2 


is preserved 

The heat-capacity curve (fig. 1) is smooth in 
range of temperature in which run 4 was made 
polymer Is considered to be completely amorph 
above 285° K. 

The remainder of the experiments were made aft 
cooling the polymer slowly to obtain a well-cryst 
lized and annealed sample. The results of 
measurements In runs 5, 6, and 7, which were mad 


on successive days, showed that, in general, afte; 
annealing the polymer, the heat-capacity val 
between 210° and 285° K became higher (see ta! 
2, and 3,and fig. 1). This indicates that the amount 
crystallization that took place during these heat 
capacity measurements was lower than in th 


measurements when the polymer was shock-cooled 
In these experiments, the downward temperat 
drifts below 200° K are considered to arise from th 
slow transition of the polymer molecules from lower 
to higher-energy states, these states being associated 
in some way with glass transformation. The dow: 
ward temperature drifts above this temperature a1 
attributed to slow melting of the crystallites. As 
relatively low heat-capacity results were obtained i 
ved heat capacities of 41 tadue - | the temperature interval 200° to 230° K, the polym 
Ryrees (ory ep rs was held within this temperature interval for 1 week 
to increase the degree of crystallinity. Runs 8 and 9 
carried out on successive days, are experiments wit! 
the material so treated The heat-capacity results 
in these measurements were not as low as found 
runs 2 and 6; however, a relatively low value was 
obtained at 227.60° K. The experiments in run 
were made to determine whether a longer cond 
tioning period would eliminate completely the low 
heat-capacity values. In spite of the 3-weeks 
conditioning, the experiments of run 12 yielded low 
values. In this run,the polymer was probably not as 
well crystallized as expected. The interpretation Ol 
the observed results with this polymer is mad 
difficult by the closeness of the crystallization to the 
glass-transformation temperature. In the crystal- 
lization interval, probably the melting and crystal 
lization processes occur simultaneously. The diree- 
tion of the temperature drift is an indication of the 
predominating process at the time of the measure- 
ment. Also, the direction of the drift is dependent 
as in the glass-transformation temperature rang 
upon the previous heat treatment, the rate at whi 
the new condition was reached from the former, at 
the time lapse between the attainment of the n 
condition and the experimental observation. 
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The annealing process, which increased the erystal- 
nity in the polymer, yielded lower heat-capacity 
ilues in the temperature range below 200° K. The 
ference between the results of the shock-cooled 
d annealed polymer is about 1 percent at 185° K 
d about 0.5 percent at 80° K. The heat-capacity 
ilues given In table 4 were calculated from the 
sults of the annealed polymer. The values at 15° K 
d below were obtained by extrapolation, using the 
lebve function, 


C=0.2489D ( _ ), 


hich was fitted to experimental values between 17 

d 25° K. The heat-capacity values in the range 
10° to 295° K were estimated from a heat-capacity 

mperature curve, in which the values in the crystal- 
ation range were taken along the broken curve 
:ypothetical heat-capacity curve for amorphous 
»polymer), as shown in figure 1. The estimated 
ilues are given in table 4. The broken curve was 


btained by linear extrapolation of the experimental 
alues in the temperature range 290° to 330 kK. 


bove the glass-transformation temperature, the 


eat-capacity curves for Hycar 0.R.—15 [1] and GR 


’} were found to be closely linear. The results with 


he GL-658 copolymer, which was found to give no 


rystallization, are also closely linear between 250° 
nd 330 
transformation temperature, however, would result 


K: extrapolation below 250° to the glass- 
| 


a deviation of 4 percent at 200° K 


BLE 4 Heat capacity, enthalpy, and entropy of 41° | 
hutadiene (90)-styrene (10) copolymer X—45 
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The enthalpy values given in table 4, except in the 
range 190° to 295° K, were obtained by evaluating 
the expression 

“7 
H,—H CdT l 


Hy, is the enthalpy at absolute zero of temperature 
The enthalpy change over the temperature interval 
190° to 295° K was computed from the experimental 
heat-capacities obtained in run 7 by summing the 
product of individual net heat capacities and the 
corresponding temperature rise. At other tempera- 
tures, except below 15° K, the expression in eq (1 
was evaluated by numerical integration, using four- 
point Lagrangian integration coefficients [11]. Be 
low 15° K, the enthalpy was evaluated by using the 
Debye function given above 

As previously stated, in the crystallization interval 
the hypothetical amorphous copolymer would proba- 
bly follow the broken heat-capacity curve By 
integrating along this curve (see fig. 1 and table 4 
the change in enthalpy of the “amorphous” copoly- 
mer between 190° and 295° K was found to be about 
180 abs j g Comparing this value with the 196.46 
and 198.97 abs j g~' obtained in runs 7 and 9, 
respectively, for the same temperature interval, the 
heat of fusion per gram would be about 16 to 19 abs j 
The range of these values is close to 16.71 abs Jz 
obtained by Bekkedahl and Matheson [10] with 
natural rubber. As shown by the difference in the 
values obtained in runs 7 and 9, these values are 
dependent upon the heat treatment. As the copoly- 
mer was annealed for a longer period prior to run 
9 than run 7, the larger value in run 9 seems con- 
sistent 

Entropy values given in table 4 were obtained in 
a similar manner as were enthalpies by evaluating 
the thermodynamic expression 


*TCdT 
T . — 


So is the entropy at absolute zero of temperature 
Between 190° and 295° K, the entropy change was 
computed using the experimental data obtained in 
run 7 by summing the individual CA7)/7 values 
C is the net heat capacity, and AJ’ and T are the 
corresponding temperature increase and midtem- 
perature, respectively 

The enthalpy and entropy values given in table 4 
were evaluated from the results of the experiments 
with the annealed polymer. No attempt was made 
to make a similar table of the results of shock 
cooled experiments 


4.2. GL-658, 122° F Copolymer 


The heat. treatments and the subsequent heat- 
capacity experiments with the GL-658 copolymer 
were carried out in a similar manner as with the 
X—454 copolymer. The details regarding the heat 


treatments, temperature range of the measurements, 
and the drift observations with this copolymer are 








Ternperature, * 


summarized in table 5 The observed heat-capac ity 


values are given in table 6 and plotted in hgure 2 


I made afte! 


The experiments in runs | and 3 were 
rapidly cooling the copolyme! In these Measure- 
ments upward temperature drifts were observed in 


the temperature range from about 167° to 195° K 

and no drift above 195° K. The results of runs 1, 2, 

and 3 show that the glass transformation occurs at 

about 190° to 195 

this copolymer The 
] 


K and that crystallization does 
not occur mn heat-capacity 
results (run 4 
slowly, showed downward temperature drifts from 
189° to 195° K and no drift above 195° K The 


downward drift in temperature is considered to arise 


from slow changes in the polymer to higher states of 


experiments the heat-capacity 
interval 


energy In these 
curve at the glass-transformation 
few «lk vyrees soonel than the curve obtained when the 
copolymer was rapidly cooled The glass-trans- 
temperature is taken to be 193° K, the 


rises a 


formation 
temperature at the midpoint of the upward sweep in 
the curve Between 90° and 165° K, the heats 
capacity re sults seem essentially the same, regardles- 
of the heat treatment. The maximum spread in the 
over-all heat-capacity results is 0.2 percent at the 
higher temperatures (above 195° K From 165° K 
to the glass-transformation temperature the results 
show greater scattering, as mu h as ] percent The 
heat-capac itv values show no clear correlation with 
heat treatment because of the obscuring effect the 
temperature drifts have on the heat capacity At 
K the temperature was drifting upward after 
Thus, no further attempt was made in 
reach thermal equilibrium 
The temperature measurements were started at 
times (about 6 to 7 min after the end of heating 
usually required for thermal equilibrium with simpler 
Experience with the X-454 copoly- 
recently investigated 


177 
x) min 


other experiments to 


materials [5,8] 
mer and the other polymers, 


at the Bureau, shows that the waiting period for 
long. The 


thermal equilibrium can be extremely 


runs 5, 6, and 7) below 90° K 


heat-capacity results 
do not show any unusual effects 


obtained after cooling the copolymer 
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The results of the heat-capacity measurements 
were used to obtain heat-capacity values at 5-deg 
intervals from 0° to 330° K (table 7). The values 
at 15° K and below were obtained by extrapolation 
using the Debye function 


SS.7 


C 0.2633 D ( 7 


which was fitted to the experimental results obtained 
in the temperature range 17° to 25° K. 
The effect of cooling rate upon the change in en 
| thalpy between 175° and 210° K was computed fron 
| the heat-capacity data. The experiments (runs 
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and 3) made with the material cooled rapidly yielded 
17.61 and 47.75 abs j g“', respectively, giving an aver- 
age of 47.68 abs j g The result with the polymer 
slowly cooled (run 4) was 48.98 abs j g These re- 
sults show that when the polymer is slowly cooled, 
the enthalpy change, between the temperature limits 
viven, Is greater by about 1.30 abs } g 3 

The enthalpy values (table 7) at 5-deg intervals 
from 0° to 330° K were obtained in a similar manner 
as for the X-454 copolymer by evaluating eq 1. For 
the temperature interval 175° to 210° K, the average 
enthalpy change for the three experiments (1, 3, and 
4) was used in constructing the table. No attempt 
was made to evaluate the enthalpies at other temper 
atures within this interval 

The entropy values were evaluated in a similar 
manner as were those of X-454 copolymer. These 
are given in table 7. The entropy change over the 
temperature interval 175° to 210° K used in the com- 
putations is the average of the three experiments 
previously mentioned 


4.3. Reliability of Results 


The reliability of the heat-capacity measurements 
with the copolymers is difficult to evaluate because of 
the nonreproducibility of the physical state. Camp- 
bell and Allen [12] have shown in photomicrographs 
similarity in spherulite formation when the polymer 
was not raised too far above its melting temperature 
and recrystallized. The scattering in the heat- 
capacity results of X-454 copolymer does not seem 


to support this, however. The heat-capacity meas- 
urements with more normal materials [5, 8] in a sim- 
ilar calorimeter yielded results considered to have an 
error of 0.2 percent The results with the X-454 
copolymer show the heat capacities of the two series 
of experiments to differ as much as 1.2 percent just 
below the glass-transformation temperature. With- 
in the same series of experiments, in which the ma- 
terial has been treated very similarly, the results are 
within 0.2 percent, except near the transformation 
temperature. The results above the crystallization 
temperature (that is, above 285° K) in the amorphous 
region are believed to have the error of 0.2 percent 
Below the glass-transformation temperature to 50° K 
the error of the results in the same series of experi- 
Below 50° K. 
the error increases to several percent at 17 K 

The results with the GL-658 copolymer show as 
much as l-percent scattering in the temperature 
range 165° K to the glass-transformation tempera- 
ture Between 90° and 165° K the results are essen- 
tially the same within +0.2 percent, regardless of the 
heat treatment Above the 
temperature at which the heat-capacity results lie on 
a smooth curve, the results are considered to have an 
error of 0.2 percent Between 50° and 165° K, the 
error is believed to be 0.3 percent and between 165° K 
and the glass-transformation temperature, to be 0.5 
to 1 percent Below 50° K the error increases to sev- 
eral percent similar to the X-454 copolymer 

Except for the transition intervals, the tabulated 
tables 4 and 7) enthalpy and entropy for the two 
copolymers were obtained, as previously shown, from 
the smoothed heat capacities These thermal quan- 
tities are dependent upon the accuracy of the heat 
capacity from which they were evaluated. Consider- 
ing the various sources of inaccuracy, the error in the 
tabulated values of these quantities is believed to be 
0.5 percent 


ments is considered to be 0.3 percent 


glass-transformation 


5. Discussion 


The results of these experiments show that the in- 
crease in the polymerization temperature from 41° to 
122° F eliminated crystallization effects in the 90/10 
butadiene-styrene copolymer. This behavior, in 
which the crystallizabilitv of polymers is increased 
with decreasing polymerization temperatures, has 
been observed by other investigators. Beu, et al 
[13] and Beu [14], using an X-ray method, found such 
a behavior in the crystallinity of stretched polybuta- 
dienes and butadiene-styrene copolymers prepared at 
temperatures from —20° to +55° C. Lueas, et a 
[15], using a dilatometer method, studied relative 
erystallizability and molecular regularity in polybuta- 
dienes, polychloroprenes, and  butadiene-styrene 
copolymers prepared at temperatures ranging from 

33° to +50° C, and found a similar behavior in the 
crystallizability of these polymers. More recently 
Campbell and Allen [12] used a polarizing microscope 


‘ For these experiments a true probat ror 

he values given are estimate che 

curacy from all known sources The authors estin 
chance that the error is no larger than that 
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to show greater crystallizability in polybutadienes | expansion measurements is closely associated w 


prepared at lower temperatures. They found 41° F | the brittle-point temperature, below which rub 
90/10 butadiene-styrene copolymer to show a small | polymers lose rubber-like elasticity. Simila: 
ervstallinity dielectric-relaxation time [19] is shown to iner 


3ekkedahl and Matheson [10] measured the heat | steeply as temperature is lowered through the gla 
capacity of natural rubber and reported a fairly | transformation temperature. There are many s 





sharp melting point, giving the melting temperature | gestions [20, 21, 22] as to how molecular structur: 
as 284° K. Wood and Bekkedahl [16] found the | the polymer affects the glass-transformation temp: 
crystallization temperature of unvulcanized natural | ture. The size of the side groups, symmetry of 


rubber to extend from 220° to 290° K. Their polymeric units, degree of cross-linking, and eas¢ 
measurements also show the melting temperature to | rotation about the carbon-carbon bonds as might 





depend upon the crystallization temperature. The | influenced by the intramolecular structure a 
results of the measurements with the X—454 copoly- | neighboring molecules are some of the factors t] 
mer reported herein show the crystallization tem- | seem to affect the glass-transformation temperatu 
perature to extend from 210° to 285° K, with no | The transformation temperature in styrene-divin 
indication of any sharp melting temperature. There | benzene copolymers [21] has been shown to incre: 
is, however, at about 250° K a maximum (fig. 1) in | when the number of monomer units between cro 
the heat capacity. In recent unpublished heat- | links becomes less than 300. The crystallizati 
capacity investigations at the Bureau with 41° and | process would be expected to have a similar eff: 
122° F polybutadienes, the crystallization tempera- | [21], if the degree of crystallization is such that tl 
ture ranges were found to be from 200° to 295° K | number of monomer units between crystallit 
and 200° to 270° K, respective ly becomes smaller than the above figure The hea 
In the heat-capacity experiments by Bekkedahl | capacity measurements [10] with natural rub! 
and Matheson [10] the glass-transformation tempera- | show that this transformation temperature is hig! 
ture was reported to be 199° K. Thermal expansion | when the material is crystallized, although attentio 
measurements, in which the coefficient of thermal | has not been called previously to this fact. Th: 


expansion rises steeply with temperature at the | results with the X-454 copolymer show a similar ris 
vlass-transformation temperature in a similar manner | in the transformation temperature when the samp 


as the heat-capacity curve, show glass transformation | was annealed \s previously pointed out, tl 
in natural rubber to occur at 200° K [17] Ferry and results are, however, dependent somewhat upon th: 
Parks [18] investigated the heat capacity and thermal | time allowed before temperature measurement. O 
expansion of polyisobutylene and found the glass- | the other hand, the experiments with the 122° | 
transformation temperature in both measurements | copolymer, which does not crystallize, seem to shi 

to be 197° K. As pointed out by Kauzmann [19], | that the steep increase in the heat-capacity cur 
this indicates a ¢ lose alliane e of motions in polymers occurs a few degrees lower in the slowly cooler 
involving energy and volume. Also, glass trans- | experiments. In comparing the glass-transformatio: 


formation as observed in heat capacity and thermal- temperatures the experimental results with the two 
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Figure 3 Comparison of the heat-capacily curves of various 
polymers. 






Natural rubber [10); 2, GR-S [2]; 3, Hycar O. R.-15 [1]; 4, polyisobutylene 
) ; F 90/10 butadiene-st yrene copoly mer; 6, 41° F polybutadiene; 7, 122° I 


y butadiene, and 8, 122° F 90/10 butadiene-styrene copolymer 
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ed w volymers indicate that the change in the poly- 
rub rization temperature from 122° F (50° C) to 41 
milar 5° C) raised the transformation temperature from 
iInere to 200° K. This behavior is similar to that 
e gla ind with the recently investigated 122° and 41° F 
ny s vybutadienes, in which the glass-transformation 
clure nperatures obtained are 187° and 195° K, re- 
‘Mm pe tively These results seem accordant on the 
y of t sis that a higher polymerization temperature 
east lds a polymer of lower regularity [12, 13, 14, 15] 
ight a more open intermolecular structure to permit 
re a ater ease in the rotation of the polymer segments 
rs tl \s previously mentioned, crystallization in 
ratu ases the glass-transformation temperature, and 
divin lowering of the polymerization temperature 
here: reases the ery stallizability These effects are all 
1 Cro sistent with the experimental results 
izati Puckett [22] and Uberreiter [23] have shown that 
r eff styrene-butadiene copolymers the styrene hinders 
rat ti rotation, and that the increase in the bound- 
tallit rem content raises the glass-transformation 
e heat iperature. The glass-transformation temperature 
rubb and 187° K for the 41° and 122° F polybuta- 
hig! nes, 200° and 193° K for the 41° and 122° F 
tentiol tadiene-styrene (8.58 percent) copolymers re- 
Th rted here, and 237° K very recently obtained at 
lar Sureau for the 122° F b itadiene-styrene (42.98 
samp cent) copolymer are in conformity with the 
t, the sults of the above investigators 
on thi The heat-capacity curves for polyisobutylene [18], 
t 0 tural rubber [10], Hyear O. R.-15 [1], GR-S [2], 
22° | and 122° F polybutadienes, and 41° and 122° F 
» sho 10 butadiene-styrene copolymers are compared in 
cury vure 3 It is interesting to note that the heat 
cooled pacities per unit mass for these polymers lie within 


narrow band The heat of the two 
opoly mers reported merge Into a single curve below 
he glass-transformation temperature; percentage 
the curves are not quite the same as 
iown by tables 4 and 7. The experimental values 
the transformation temperature for the two 
polymers differ by about 0.1 percent at 295° K, 


natiol 
ne two 


capacities 


howey el 


ove 


creasing to about 0.7 percent at 330° K. The two 
olybutadienes, previously mentioned, show similar 


ehavior. The glass-transformation temperatures 


or 


the various polymers are compared in the figure 
highest 


howing Hycar 0. R15 [1] to have the 


transformation temperature. 
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